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ABSTRACT 

We present a near-infrared photometric search for unresolved substellar companions and de- 
bris disks around white dwarfs in the UKIRT Infrared Deep Sky Survey Large Area Survey. 
We cross-correlate the SDSS DR4 and McCook & Sion catalogues of white dwarfs with 
the UKIDSS DR8 producing 3109 and 163 unique matches respectively. Cooling models are 
fitted to the optical photometry of a subsample of DA white dwarfs and extended to the near- 
infrared. A comparison is then made with the observed photometry to identify those stars 
with a near-infrared excess consistent with the presence of a cool companion or debris disk. 
Where present, we have estimated the approximate spectral type of any putative companion, 
or an upper limit on the temperature of a debris disk. In total we identify 14-16 new candidate 
white dwarf + very low mass stellar systems, 9-11 candidate white dwarf + brown dwarf 
systems, and 3 candidate white dwarf + debris disks. We place lower limits on the unresolved 
(< 2") companions to all DA white dwarfs and thus assess the sensitivity of UKIDSS to such 
objects. We use this result to estimate unresolved binary fractions of /wD+dL ^ 0.4 ± 0.3%, 
ZwD+dT > 0.2% and / W d+bd > 0.5 ± 0.3%. 

Key words: Stars: white dwarfs, low-mass, brown dwarfs, infrared: stars 



> 1 INTRODUCTION 



Near-infrared (NIR) surveys of white dwarfs enable the detec- 
tion and study of late stellar and substellar companions, and cir- 
cumstellar dust disks. A typical white dwarf is 10 3 — 10 4 times 
fainter than its main sequence progenitor, significantly reducing 
the brightness contrast problem when searching for cool, low mass 
secondaries and dust. In addition, the spectral energy distribu- 
tions of most white dwarfs (blue) and their low mass compan- 
ions (red) are markedly different, facilitating easy separation of 
the components in broadband photome try and ena bling straightfor- 
ward spectroscopic follow-up (e.g. iDobbie et alj|2005f) . Previous 
NIR p h otometric surveys of wh i te dw a rfs include those by 



1983) 



2000) 



Probst 



Zuckerman & Becklinl dl987 ). [Green. Ali & Napiwotzki 
Farihi, Beck lin & Zuckerman < l2005h . and the analysis 



of white dwarfs in the Two-Micron A ll Sky S urvey (2MASS; 
ISkrutskie etail 120061) by IWachter et al.1 fe003l). IWellhouse etafl 
d2005l).lHolberg & Magargall d2005l) . lTremblav & Bergeronl j2007h 
and lHoard et alj d2007l) . 

Many white dwarfs are known to have low mass M dwarf 
companions, and analysis of their population statistics allows the 
investigation of binary formation and evolution. In particular, con- 
sideration of stellar evolution suggests the possible existence of two 
distinct populations: close systems in which the secondary has sur- 



vived a phase of common envelope evolution and which may even- 
tually lead to the formation of a cataclysmic variable (CV), and 
wide pairs where the secondary has migrated out wards in respons e 
to mass-loss from the white dwarf's progenitor dFarihi et alj 2006). 
One particularly interesting case are the magnetic white dwarfs. 
Although hundreds of magnetic CVs (Polars) are known, no con- 
vincing progenitor system has been identified dLiebert et al.ll2005h 
despite exhaustive anal ysis of the 2MASS dataset for these objects 
dWellhouse et alj|2005h . 

Searches for substellar companions to white dwarfs addi- 
tionally allow the investigation of the known deficit of brown 
dwarf companions to main sequence stars (< 10AU: 0.5%; 
iMarcv & Butieil 2000) The closest brown dwarf + white dwarf 
binaries might also represent either another channel for CV evo- 
lution or the end state of CV evolution, in which the secondary 
has become highly evolved through mass transfer. In close de- 
tached binaries, the brown dwarf is expected to be irradiated by 
the white dwarf's high UV flux, possibly leading to substan- 
tial temperature differences between the "day" and "night" hemi- 
spheres. Such systems can provide laboratories for testing mod- 
els of irradiated "h ot Jupiter" atmospheres (e.g. HD 189733b; 
Knuts on et al. 1 120071) . Howe ver, detached brown dwarf compan- 
ions to white dwarfs are rare. lFarihi. Becklin & Zuckerman! ( 2005) 
suggest a binary fraction of /wD+dL < 0.5% for L dwarfs 
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(We re-assess this limit later in this work). Radial velocity and 
proper motion surveys, and searches for NIR excesses have so 
far found o nly five spectroscopically co nfirmed examples: GD 165 
(DA+L4. iBecklin & Zuckermanl 1 19881). GD 14 00 (DA+L6 - 7; 
iFarihi & Christopher! |2004|; iDobbie et alj|2005h WD 0137 - 349 
(DA+L8; iMaxted et al.1 l200d iBurleigh et all 120061) . PHL5038 
(DA+L8; ISteele et allbOOgh and LSPM 1459+0857 (DA+T4.5; 
iDav- Jones et alj201lh . GD 165, PHL5038 and LSPM 1459+0857 
are wide orbit binaries with projected separations of 120 AU, 
55 AU and 16,500-26,500 AU respectively, whereas GD 1400 and 
WD 0173— 349 are in much closer orbits with periods of 10 hours 
(Burleigh et al. 201 1, in prep) and 1 16 mins respectively. 

Infra-red observations have recently revealed the pres- 
ence of metal-rich debri s dust disks around a small num- 
ber of white dwarfs (e.g. Becklin et al. 2005; Kilic et all 



2005; 



Kilic et alj 120061; iKilic & Redfield [20071: Ivon Hippel etai] |2007| ; 
Jura. Farihi & Zuckermanl 120071) . Indeed, the large NIR ex- 
cess in G29— 38 w as first identified twenty years ago by 
IZuckerman & Becklinl dl987l) . and Spitzer mid-infrared spec- 
troscopy ha s now revealed tha t this dust disk is composed largely 
of silicates ( Reach et al. 2005). The favoured explanation for the 
origin of this material is the tidal disruption of an aste roid that ha s 
wandered into the Roche radius of the white dwarf Jjur3l2003h . 
With one notable exception, these disks have been found around 
relatively cool white dwarfs whose atmospheres are polluted with 
heavy elements. These metals should sink from the photospheres 
on a timescale of days, so the presence of close, orbiting dust disks 
provides a reservoir for on-going low level accretion. The excep- 
tion is the central star of the Helix nebula. This hot, extremely 
young whi te dwarf has a mid-IR excess emission identified with 
Spitzer bv lSu et aT I ( 12007b . which might be explained by the pres- 
ence of a ring of colliding cometary material at an orbital distance 
of ~ 30 AU. Therefore, dust disks around white dwarfs may have 
a variety of origins linked to the post-main sequence evolution of 
both the central stars and surrounding planetary systems, and large 
IR surveys will help establish how common they are. 

In this paper we present the first results from studying white 
dwarfs detected in the eigth data release (DR8) of the UKIRT In- 
frared Deep Sky Survey, the largest and deepest survey of the near- 
infrared sky yet undertaken. 



2 THE UKIDSS SURVEYS 

The UKIRT Infrared Deep Sky Survey (UKIDSS ) is the NIR coun- 
terpart of the Sloan Digital Sky Survey (SDSS; lYork et al.ll2000h 
and is several magnitude s deeper than the previous 2MASS sur- 
vey dSkrutskie et alj|2006l) , although unlike the latter UKIDSS will 
cover only ~ 20% of the sky. UKIDSS actually consists of a set 
of five complementary wide/shallow, narrow/deep (extra-galactic), 
Galactic cluster and Galactic Plane surveys. A full description of 
these surveys, the science goals, observations a nd data processing 
and products are given in lLawrence et alj 00071) . Briefly, UKIDSS 
uses the 0.21 deg 2 field of view Wide Field Camera (WFCAM; 
ICasali et al]|2007l) on the 3.8m UKIRT telescope on Mauna Kea, 
Hawaii. The UKIDSS surveys began in May 2005 and are expected 
to take 7 yr to complete. The survey of most relevance to this work 
is the Large Area Survey (LAS; PI: S. Warren), which aims to cover 
~ 4000 deg 2 of the Northern Sky coincident with th e SPSS. The 
LAS makes observations in the YJHK filters (see iHewett et al.1 
2006 for a description of the UKIDSS photometric system) to a 5a 
depth for point sources of Y « 20.2, J « 19.6, H « 18.8 and 



K x 18.2. At K the LAS is 2.7 magnitudes deeper than 2MASS. 
UKIDSS data are being released to the community in stages and 
are public immediately to astronomers in all ESO states, before 
being made available to the rest of the world after 18 months. Uni- 
formly processed images and catalogues are available through the 
WFCAM Science Archiv^B Following an the r elease of scienc e 
verification data and an Early Data Release (E DRuDve et al j 2006). 
UKIDSS made the First Data Release (DRL lWarren et al.ll2007bh 
in July 2006 and the Second Data Release (DR2; Iwarren et all 
12007 at) in March 2007. A further 6 data releases haven subsequently 
been released (DR3-DR8) approximately once every 6 months. 
The most recent release, DR8, is inclusive of all previous data re- 
leases, and covers 2445 deg 2 in all four YJHK filters, represent- 
ing ~64% of the final survey. 

The LAS observations are carried out in the YJHK filters 
with the goal of reaching a uniform depth in all fields. The mini- 
mum schedulable blocks (MSBs) for the LAS use a pair of filters; 
either Y J or HK. If an area has been covered by one pair an ef- 
fort is made to follow up with the remaining pair within one month. 
However, conditions do not always allow for this and observations 
depend on a number of constraints (Dye et al. 2006). Hence, before 
completion a small number of observations in each data release will 
only have measurements in either the Y J or HK bands only. It is 
also entirely possible that any single measurement will not pass 
data quality tests and be omitted from the database. Therefore, a 
number of objects may have a selection of Y, J, H or K magni- 
tudes, but not all four. 



3 IDENTIFICATION OF WHITE DWARFS IN THE 
UKIDSS LAS DR8 

We have cross-correlated the UKIDSS LAS DR8 with two cata- 
logues of white dwarfs: the on-line August 2006 version of the 
McCook & Sion catalogue (hereafter MS9 9) of Spectr oscopically 
Identified White DwarfsQ (see iMcCook & Sionlll999l for a fuller 
description) which contains 5557 stars, and the 9316 spectroscopi- 
cally confirmed white dwarfs catalogued by Eisenstein et al. ( 2006) 
(hereafter EIS06) from the SDSS Data Release 4. The MS99 cata- 
logue includes the brightest white dwarfs known. The white dwarfs 
identified in the SDSS are all fairly hot (T e g > 8000 - 9000 K), 
and are constrained by the upper brightness limit (</ w 15) and 
lower sensitivity limit (</ ~ 21) of that survey. 

Inevitably, there is considerable overlap between the cata- 
logues. However, we note that the majority of SDSS white dwarfs 
are too faint to have be en detected in the 2MASS survey (for exam- 
ple, |HoMdetal]|2003 excluded all objects J < 15.8, H < 15.1, 
K B ^ 14.3). Therefore, UKIDSS allows us for the first time to 
investigate the NIR properties of the large number of new white 
dwarfs discovered in the SDSS, as well as significantly reducing the 
measurement errors on the NIR photometry of the brighter white 
dwarfs detected by 2MASS. 

Many of the relatively bright white dwarfs have large and un- 
measured proper motions, meaning their astrometry is often unreli- 
able. For the MS99 white dwarfs, we have made use of a set of im- 
proved co-ordinates kindly supplied by Jay Holberg (private com- 
munication) in our cross-correlation with UKIDSS DR8. The SDSS 
stars, being generally fainter and more distant, have small proper 
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Table 1. Number of each white dwarf spectral type detected in the cross 
correlation of UKIDSS LAS DR8. 



Spectral Type 


No. 


Spectral Type 


No. 


Spectral Type 


No. 


DA 


2306 


DO 


10 


DBAZ 


1 


DA+M 


350 


DBZA 


3 


DAO 


1 


DC 


94 


PG1159 


5 


DBZ 


1 


DB 


153 


DZA 


6 


DO+M 


2 


DQ 


31 


DBH 


3 


DAZ 


2 


DZ 


37 


DBZ 


4 


DZ+M 


1 


DAH 


19 


DAB 


5 


DBA+M 


1 


DBA 


37 


DH 


2 


DAH+DA 


1 


WD+M 


13 


DAOB 


2 


DAH+M 


1 


DB+M 


6 


DACV 


2 


DBO 


1 


DC+M 


5 


DZA+M 


1 


DQ+M 


1 










DZB 


1 



motions and generally reliable astrometry. We require all matches 
have an S/N ^ 5, corresponding to J ^ 19.56, H ^ 18.81 
and K ^ 18.23 dWarren et alj|2007ah . We performed the cross- 
correlation using the on-line "CrossID" tool at the WFCAM Sci- 
ence Archive. We used a pairing radius of 2", and took the near- 
est NIR source to the white dwarf's co-ordinates as being the best 
match. In general, the UKIDSS and white dwarf co-ordinates agree 
to better than 2" and we are confident that the correct match has 
been made. A few SDSS "white dwarfs" had very red UKIDSS 
colours: an inspection of the corresponding Sloan spectra revealed 
these to be quasar contaminants of the EIS06 catalogue. 

The cross correlation produced 3109 matches with the EIS06 
catalogue, and 163 exclusive matches with the MS99 white dwarf 
catalogue. The corresponding spectral types for the white dwarfs 
(as assigned by EIS06 for the SDSS DR4 white dwarfs, and from 
literature for the MS99 white dwarfs) are listed in TableQ] It should 
be noted that the EIS06 catalogue is not a complete sample of white 
dwarfs, but the result of the placement of SDSS spectroscopic fibers 
according to complex and largely non-stellar considerations. Like- 
wise, the MS99 catalogue is a heterogeneous sample. Thus, a num- 
ber of white dwarfs outside of this sample may still be present in 
the area covered by the UKIDSS surveys. 



4 ANALYSIS 

4.1 The Near-Infrared Two Colour Diagram 

In order to make an immediate and rapid search for white dwarfs 
with a NIR excess, indicative of a low mass companion or dust disk, 
the UKIDSS DR8 white dwarfs were plotted on the (J - H,H - K) 
two-col our diagram shown i n Figure Q] This method was first ap- 
plied bv lWachter etail d2003h in their analysis of white dwarfs de- 
tected by 2MASS. Then, b ased on Wachter et al.'s work and the- 
oretical colour simulations, IWellhouse et"ai] (2005) further devel- 
oped the technique by splitting the two-colour diagram into four 
distinct regions; (I) single white dwarfs, (II) white dwarfs with 
late-type main-sequence companions, (III) white dwarfs with lower 
mass, possible brown dwarf companions and (IV) white dwarfs 
contaminated with circumstellar material. Wellhouse et al.'s regions 
have been plotted in th e two-colour diagram in Figure[T] 

Hoard ct al] d2007h argue this colour selection method allows 
a simple, fast and efficient means to exploit large, homogeneous 
survey data to provide an initial list of candidate stars with an inter- 
esting NIR excess. Hoard et al. (2007) confirmed 80% of the white 



Figure 1. J— H v H— K colour-colour diagram for the EIS06 white dwarfs 
(top) and the MS99 (bottom) white dwarfs with complete JHK photome- 
try present in the UKIDSS archive. The error bars have been omitted from 
the EIS06 diagram for the sake of clarity. 
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dwarf + main-sequence binaries reported by Wachter et al. (2003) 
with a false positive rate (i.e. stars rejected as misidentification by 
Wachter et al. 2003) of 20%. However, they failed to identify any 
new candidate brown dwarf companions or dust disks. This is prob- 
ably in part due to the excess emission from additional unresolved 
sources being small compared to the errors on the 2MASS data. 
However, it could also be due to some white dwarfs with compan- 
ion brown dwarfs or dust disks occupying regions (I) and (II) of 
the two-colour diagram, making them indistinguishable from sin- 
gle white dwarfs or the large population of white dwarfs with main 
sequence K and M dwarf companions. 

To investigate this further, white dwarfs with known brown 
dwarf companions and dust disks have been plotted in a NIR two- 
colour diagram in Figure [2] using photometry from 2MASS and 
other sources in literature (See Tables [2] and 0. Clearly only eight 
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Figure 2. NIR colour-colour diagram of white dwarfs with confirmed low- 
mass companions (circles) and circumstellar disks (sq uares). The data are 
given in Tables|2]and[3] The four regions denned by Wachter et all fe003l) 
and lWellhouse et aU 120051) and described in the text are shown for compar- 
ison with Figure 1 . 



of the fourteen white dwarfs plotted lie within regions where fur- 
ther investigation would have been solicited using Wellhouse et al.'s 
(2003) regions. Indeed only a single dust disk (GD56) and a single 
white dwarf brown dwarf binary (SDSS J121209.31+013627.7) lie 
within the regions expected of such systems. In this case we can 
crudely conclude that the colour selection method is likely to eas- 
ily distinguish ~ 50% of such systems. However, more examples 
would be required to accurately constrain such statistics. 

It should also be noted that this method can only be used to 
analyse white dwarfs which have complete JHK photometry. In 
Section|2]we noted that some stars will only have either Y J or HK 
photometry at this time in UKIDSS and so can not be included in 
this analysis. A massive 70% of the white dwarfs in this sample are 
missing one of the required J, H or A-band magnitudes (Table|4j 
needed for this type of analysis. 

Other limitations of the co l our s election method are dis- 
cussed by iTremblav & Bergeron (2007]), who critically analysed 
the 2MASS photometric data for white dwarfs and the analyses of 
Wachter et al. (2003) and Wellhouse et al. (2005). They concluded 
that the method is appropriate but should be interpreted with cau- 
tion near the detection threshold where larger errors are to be ex- 
pected. They demonstrated that there exists a number of single stars 
with large photometric uncertainties that contaminate regions (III) 
and (IV) of the colour-colour diagram. Indeed, It can be seen in Fig- 
ure [TJ that there are a number of stars with errors bars that extend 
from regions (III) and (IV) back into regions (I) or (II). 

Tremblay & Bergeron (2007) demonstrate that a more thor- 
ough technique is to compare observed NIR photometric fluxes 
with those predicted by m odel atmospheres (as also adopted by 
Holberg & Magargal 2005). This method first requires accurate ef- 
fective temperatures and surface gravities, which are determined 
from high S/N optical spectroscopy. The model spectral energy dis- 
tribution (SED) is then normalised to the observed V or g' magni- 
tude of the white dwarf. Predicted photometry can then be calcu- 
lated by folding the synthetic SED through the appropriate filter 
transmission profiles and a direct comparison made. However, this 



Table 4. The combinations of UKIDSS data available for the entire white 
dwarf sample 



Bands Available 


No. 


Bands Available 


No. 


YJHK 


926 


JH 


28 


YJH 


477 


JK 


1 


YHK 


28 


HK 


67 


YJK 


43 


Y 


708 


JHK 


56 


J 


81 


YJ 


790 


H 


29 


YH 


22 


K 


4 


JK 


2 







technique is limited to those stars with known effective tempera- 
tures and surface gravities. Many of the MS99 white dwarfs do have 
accurately determined temperatures and gravities. The SDSS DR4 
white dwarfs do have published atmospheric parameters (EIS06), 
however these are determined in an automated fashion and users of 
the catalogue are warned to be cautious when adopting these val- 
ues. 

The two colour diagram of the UKIDSS white dwarfs (Fig- 
ure [T} clearly separates out those stars with optical companions 
into the bottom of region (II). The majority of stars are in region 
(I) potentially indicating that these are single white dwarfs. There 
are then a small number of stars in regions (III) and (IV) which 
would warrant further investigation. It should be noted that a num- 
ber of the clear outliers (e.g. the MS99 object near the top edge 
of the diagram) are classified in the UKIDSS database as merged 
sources. This could indicate contamination, perhaps by an back- 
ground galaxy or star. It can also indicate that a companion is be- 
coming resolved in the NIR, as is most likely the case for e.g. 
SDSS 1020538. 12+005835.3 which clearly has a companion in the 
optical (Figure [3} and UKIDSS classes as a merged source, as can 
be seen from the science images (Figure[4]l. The UKIDSS merged- 
Class statistic (-1 for a point source, +1 for an extended source) 
was checked for each candidate object to assess these possibilities 
(TableUl. 

It is clear there are significant advantages to using the mod- 
elling technique over the two colour diagram; a direct compari- 
son can be made between observed and expected photometry de- 
termined from UKIDSS filter transmission profiles, and those ob- 
jects with photometry in only one pair of filters can be analysed. 
Therefore we adopt this method for our work. 



4.2 White Dwarf Modelling 

In order to efficiently process the cross-matched catalogue of white 
dwarfs we decided to model each white dwarf independently and 
then inspect each for evidence of a NIR excess. We also decided 
that only those DA white dwarfs with both H and A'-band pho- 
tometry were to be modelled, as a substellar companion will give 
only a small or non-detectable excess at Y and J, depending on the 
temperature of the components. Also, H and A'-band photometry 
together allow candidate substellar companions and blackbody-like 
dust disk emissions to be distinguished from each other. 

We have used the online version of Pierre Bergeron's syn- 
thetic colours and evo l utionary sequences for DA white dwarffl 
jBergeron et~aH 1 19951 : iHolberg & Bergeronl 120061) . This grid is 



3 http://www.astro.umontreal.ca/ bergeron/CoolingModels/ 
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Table 2. Near infra-red magnitudes of known white dwarfs with unresolved brown dwarf companions 



Name J H K Source 

WD 0137 - 349 15.681 ± 0.055 15.887 ±0.148 15.438 ±0.188 2MASS 1 

GD1400 14.923 ±0.032 14.450 ± 0.045 14.335 ± 0.064 2MASS 1 

SDSS 1121209.31 +013627.7 17.91 ± 0.05 17.676 ± 0.091 17.372 ±0.114 Farihi et al. (2008) 

PG 1234 + 482 14.977 + 0.043 14.955 ± 0.073 14.937 ±0.106 2MASS 1 

PHL5038 2 16.75 ± 0.02 16.44 ± 0.03 16.36 ± 0.04 Steele et a. (2009) 

1 Skrutskie et al. (2006) 
2 Partially resolved in UKIDSS 



Table 3. Near infra-red magnitudes of known white dwarfs with circumstellar debris disks 



Name 



H 



Source 



G29 - 38 
GD362 
GD56 
WD 2115 - 
GD40 
PG1015 
GD133 
GD16 

SDSSJ122859.93 - 



560 



161 



104032.9 



13.132 ±0.029 

16.09 ±0.03 
15.870 ±0.061 
14.110 ±0.029 
15.885 ±0.072 
16.131 ±0.085 
14.752 ±0.039 
15.799 ± 0.076 
16.897 ±0.016 



13.075 ± 0.022 

15.99 ±0.03 
15.991 ±0.129 
13.996 ±0.055 
15.897 ±0.161 
16.120 ±0.222 
14.730 ±0.051 
15.528 ±0.101 
16.856 ± 0.023 



12.689 ± 0.029 

15.85 ± 0.06 
15.440 ± 0.175 
14.022 ± 0.061 
15.621 ± 0.214 
16.003 ± 0.216 
14.611 ± 0.105 
15.358 ± 0.163 
16.442 ± 0.038 



2MASS 1 
Becklin et al. (2003) 
2MASS 1 
2MASS 1 
2MASS 1 
2MASS 1 
2MASS 1 
2MASS 1 
UKIDSS DR5 



1 Skrutskie et al. (2006) 




Figure 3. SDSS finding chart for SDSS 1020538. 12+005835.3 clearly 
showing both system components (circled). The colours have been inverted 
from the original finding chart. 




Figure 4. UKIDSS YJHK images of SDSS J020538. 12+005835.3. 



built from a loc al thermodynamic eq uilibrium (LTE) model at- 
mosphere code (Berg eron et alj 1992), which incorporates con- 
vective energy transport and hydrogen molecular opacity, up to 
an effective temperature (T c ft = 20, 000 K) where non-LTE 
effects are still negligible and the stellar atmosphere is com- 
pletely radiative. For pure hydrogen white dwarfs with T e a > 
20, 000 K, the grid then switches to using the plane-parallel, hydro- 
static, non-LTE atmosphere and spectral syn thesis codes TLUSTY 
(v200. lHubenvll988l ; lHubenv & Lanzll 19951) and SYNSPEC (v49; 
Hubeny & Lanz, ftp:/tlusty.gsfc. nasa.gov/synsplib/synspec). The 
original model grid extends from 1500 < T c ff < 100, 000 K and 
7.0 <logg < 9.0 for DA white dwarfs. This was then linearly 
extrapolated to 6.5 <log(? < 9.5 to include the lower and higher 
mass white dwarfs, but to preclude any possible sub-dwarf contam- 
ination of the sample. 

This provided theoretical absolute magnitudes for the SDSS i' 
and the 2MASS JHKs filters, as well as an estimated distance to 
the white dwarf and cooling age of the system. The 2MASS magni- 
tudes were converted to th e corresponding U KIDSS filters using the 
colour transformations of lCarpenterl d200ll) . Out of the 2677 white 
dwarfs classified by EIS06 and MS99 as DA white dwarfs, 811 had 
both H and if-band photometry and fell within the extended DA 
model grid. 

We identified candidate stars as those showing at least a > 3a 
excess in the UKIDSS i\"-band when compared to the model pre- 
dicted values. 313 such excesses were identified, with 274 of these 
accounted for as previously identified binaries from an optical ex- 
cess in SDSS (or other optical data). This left 39 candidate white 
dwarfs with an apparent NIR excess. The UKIDSS photometry for 
these is given in Table [5] with additional parameters given in Ta- 
ble© 

4.3 Identifying Sub-Stellar Companions 

For each candidate white dwarf we first examined the SDSS opti- 
cal spectrum for any previousley undetected evidence of a low mass 
main-sequence companion. An increase in flux of the spectral en- 
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Table 5. UKIDSS YJHK magnitudes of the candidate white dwarfs showing near-infrared excesses. 



SDSSJ 


WD 




Y 




J 




H 




K 


003902.47-003000.3 


_ 


18 


.70±0.06 


18. 


37±0. 


08 


17.86±0.11 


17 


,41±0.10 


003923.04+003534.7 


_ 


18 


.44±0.04 


17 


87±0. 


03 


17.48±0.05 


17 


04±0.06 


012032.27-001351.1 




19 


.25±0.08 


18 


34±0. 


07 


17.82±0.07 


17 


.42±0.08 


013532.98+144555.8 


0132+145 


16 


.49±0.01 


16 


,24±0. 


01 


15.86±0.01 


15 


.56±0.02 


032317.00-002612.7 


_ 


15 


03±0.00 


14 


63±0. 


00 


14.08±0.00 


13 


.81±0.01 


033209.25+010722.1 




17 


.45±0.02 


17 


47±0. 


02 


17.75±0.10 


17 


65±0.10 


034221.55+005345.6 


_ 


16 


.48±0.01 


16 


05±0. 


01 


15.45±0.01 


15 


.31±0.01 


075327.05+244705.5 


_ 


19 


.17±0.06 


19 


22±0. 


07 


18.65±0.18 


18 


.01±0.16 


082705.04+284402.7 


0824+288 


13 


06±0.00 


12 


,72±0. 


00 


12.14±0.00 


11 


.84±0.00 


085956.47+082607.5 


_ 


18 


.33±0.03 




_ 




17.49±0.03 


17 


09±0.08 


090759.59+053649.7 


_ 


18 


J9±0.08 


18. 


66±0. 


08 


18.62±0.16 


18 


.40±0.19 


092648.84+102828.8 


_ 


19 


.57±0.12 


19 


19±0. 


16 


18.74±0.19 


18 


29±0.19 


093821.34+342035.6 




18 


.48±0.05 


18 


08±0. 


03 


17.65±0.08 


17 


35±0.10 


100259.88+093950.0 




19 


.30±0.08 


19 


19±0. 


13 


18.55±0.18 


18 


.55±0.22 


101642.93+044317.7 


_ 




_ 




_ 




17.73±0.09 


17 


59±0.09 


103220.19+011227.0 










_ 




17.76±0.08 


17 


.15±0.07 


103448.92+005201.4 


1032+011 








_ 




18.20±0.11 


18 


03±0.14 


103736.75+013912.2 


_ 




_ 




_ 




17.25±0.06 


17 


.12±0.07 


110826.47+092721.5 


_ 








_ 




16.44±0.02 


16 


.21±0.04 


113416.09+055227.2 


_ 




_ 




_ 




18.10±0.12 


17 


.94±0.20 


115808.44-012312.9 


1155-011 


18 


26±0.03 


17 


83±0. 


03 


17.22±0.03 


16 


89±0.06 


120503.99+070031.0 


_ 


18 


.47±0.03 


18 


,31±0. 


05 


17.76±0.08 


17 


.70±0.13 


121209.31+013627.7 


_ 


18 


,04±0.04 


17. 


85±0. 


05 


17.21±0.05 


16 


21±0.03 


122625.93+002700.2 


1223+007 


19 


.11±0.08 


19 


13±0. 


14 


18.54±0.14 


18 


02±0.18 


122859.92+104033.0 


_ 


16 


.86±0.01 


16 


89±0. 


02 


16.84±0.02 


16 


.42±0.04 


124455.15+040220.6 


_ 




_ 




_ 




18.59±0.12 


17 


.70±0.13 


125044.42+154957.4 1 


1248+161 






16 


61±0. 


02 


16.14±0.02 


15 


.00±0.01 


131955.04+015259.5 


1317+021 




- 




- 




17.79±0.05 


17 


.57±0.10 


132044.68+001855.0 


1318+005 


17 


.61±0.02 


17 


62±0. 


04 


17.76±0.11 


17 


.38±0.11 


132925.21+123025.4 




17 


.07±0.01 


16 


66±0. 


01 


16.27±0.02 


15 


.98±0.02 


134154.39+005600.8 


1339+011 


18 


.30±0.03 


17 


82±0. 


04 


17.33±0.05 


16 


95±0.06 


141448.25+021257.7 


1412+024 


17 


.77±0.03 


17 


87±0. 


04 


17.76±0.08 


17 


19±0.09 


154431.47+060104.3 




17 


.96±0.03 


17 


83±0. 


04 


17.64±0.05 


17 


64±0.10 


155720.77+091624.7 




18 


.82±0.04 


18 


82±0. 


06 


19.03±0.14 


18 


.34±0.15 


162514.88+302610.8 




18 


.09±0.02 


17 


95±0. 


03 


17.66±0.07 


17 


.72±0.13 


220841.63-000514.5 




18 


.16±0.04 


17 


48±0. 


03 


17.06±0.05 


16 


.59±0.04 


222030.68-004107.9 




16 


.95±0.01 


16 


75±0. 


02 


16.44±0.03 


16 


.36±0.04 


222551.65+001637.7 




18 


.59±0.05 


18 


51±0. 


07 


18.15±0.13 


17 


60±0.13 


231133.88+144119.9 




19 


.51±0.11 


19 


17±0. 


15 


18.68±0.16 


17 


.88±0.17 


233345.97-000843.0 




19 


.36±0.08 


18 


80±0. 


10 


18.22±0.13 


17 


,67±0.13 



Identified in Section I5~2l 



ergy distribution towards longer wavelengths or spectral features 
indicative of M dwarfs (e.g. Ha emission and/or TiO absorption 
bands) were searched for by eye. Empirical models for low mass 
sub-stellar objects were then added to a white dwarf synthetic spec- 
trum and these composites were compared to UKIDSS photometry 
to obtain an approximate spectral type for the putative companion. 

The white dwarf pure-hydrogen synthetic spectra were calcu- 
lated at the effective tempe ratures and sur face gravities given in 
Table[6]using the models of lKoesteH j2008j). The synthetic spectra 
were then normalised to the Sloan i' magnitude of the correspond- 
ing white dwarf (Table [6}. The i'-band was chosen for this work 
to avoid the strongest absorption lines in the white dwarf's optical 
spectrum. 

The empirical models for the candidate secondary stars were 
constructed usin g the NIR spectra of late M dwarfs from the IRTF 
spectral library JCushing et al JI2005I . Rayner et al. in prep.), and 
the L and T dwarfs in Sandy Leggett's L and T dwarf a rchive 
dChiu et alJl2006| ; lGolimowski et alBoollKnapp etai]|2004h . 

In order to obtain an initial estimate on an approximate spec- 
tral type of a putative companion, we calculated an absolute K- 



band magnitude for each candidate system. This was achieved by 
subtracting the model predicted A'-band flux from the observed 
UKIDSS A'-band flux, and then scaling to the estimated distance 
of the white dwarf (Table[8}. This value was compared to observed 
absolute A'-band magnitudes of M, L and T dwarfs, which show a 
trend distinctive enough t o separate the bou ndary regions of these 
3 spectral types (Figure[5j Vrba et alj |2004). Due to errors in tem- 
perature, surface gravity and therefore distance, estimates of the 
companion spectral types are likely within ± 1 spectral types of the 
best fitting composite white dwarf + companion model spectrum. 

4.4 Identifying Debris Disks 

A debris disk around a white dwarf can be identified through a sub- 
stantial excess emission in the A'-band, whilst lacking a significant 
excess in the J and ff-band. Using the following method, black- 
body models were used to obtain an upper limit on the tempera- 
ture of the potential disk: a blackbody model was generated at the 
temperature of white dwarf and scaled appropriately. This model 
was then added to the white dwarf synthetic spectrum to create a 
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Table 6. Physical parameters of the candidate white dwarfs showing near-infrared excesses. Temperatures and surface gravities are from EIS06 unless otherwise 
stated. 



SDSSJ 


WD 




g' 


T cff (K) 


log 9 


003902.47-003000.3 


- 


19 


.11+0.02 


12314±268 


7.39±0.08 


003923.04+003534.7 


- 


19 


.20+0.02 


10080±81 


7.99±0.10 


012032.27-001351.1 


- 


20 


.23+0.04 


10487±194 


7.91±0.22 


013532.98+144555.8 


0132+145 


16 


96+0.01 


8083±22 


8.08+0.04 


032317.00-002612.7 


- 


20 


.72+0.07 


22193±1446 


7.10+0.19 


033209.25+010722.1 


- 


17 


.90+0.01 


16700+89 


7.89+0.02 


034221.55+005345.6 


- 


19 


.73+0.03 


23035±684 


6.97±0.08 


075327.05+244705.5 


- 


19 


.47+0.02 


13256±321 


7.80±0.07 


082705.04+284402.7 


0824+288 


14 


.26+0.00 


- 


- 


085956.47+082607.5 


- 


19 


,13±0.02 


15623±237 


7.44±0.06 


090759.59+053649.7 


- 


18 


.85+0.01 


19474±293 


7.82±0.05 


092648.84+102828.8 


- 


20 


.47+0.04 


13214±935 


8. 11 ±0.29 


093821.34+342035.6 


- 


19 


00+0.01 


8130±59 


8.02±0.14 


100259.88+093950.0 


- 


19 


.57+0.02 


22331±435 


7.87±0.06 


101642.93+044317.7 


- 


18 


.99+0.01 


27039±276 


7.43±0.04 


103220.19+011227.0 


- 


19 


.17+0.02 


8901±93 


7.71±0.18 


103448.92+005201.4 


1032+01 1 


19 


16+0.02 


9904±109 


8.13±0.15 


103736.75+013912.2 


- 


18 


.82+0.01 


16300±221 


7.73±0.05 


110826.47+092721.5 


- 


18 


.48+0.01 


11873±176 


8.97±0.05 


113416.09+055227.2 


- 


19 


.59+0.03 


24184±650 


7.85±0.08 


115808.44-012312.9 


1155-011 


19 


.63+0.03 


~13300 2 


- 


120503.99+070031.0 


- 


18 


,97±0.07 


8688±61 


8.21±0.12 


121209.31+013627.7 


- 


18 


.24+0.01 


~10000 3 


~8.0 3 


122625.93+002700.2 


1223+007 


20 


,17±0.04 


8744±146 


8.97±0.04 


122859.92+104033.0 


- 


16 


90±0.01 


22125±136 


8.22±0.02 


124455.15+040220.6 


- 


19 


.25+0.02 


67380±4504 


7.95±0.27 


125044.42+154957.4 


1248+161 


18 


.36+0.01 


~10000 4 


- 


131955.04+015259.5 


1317+021 


18 


,71±0.01 


16800 5 




132044.68+001855.0 


1318+005 


17.69±0.01 


19649±130 


8.36±0.02 


132925.21+123025.4 




17 


.70+0.01 


13537±138 


7.48±0.02 


134154.39+005600.8 


1339+011 


19 


.14+0.02 


18032±226 


7.91±0.05 


141448.25+021257.7 


1412+024 


18 


.18+0.01 


8296+41 


8.45+0.07 


154431.47+060104.3 




18 


.51+0.01 


7946±66 


8.39±0.15 


155720.77+091624.7 




19 


.11+0.01 


21990±403 


7.67±0.06 


162514.88+302610.8 




18 


.32+0.03 


74667±3097 


7.88±0.13 


220841.63-000514.5 




19 


.19+0.03 


9224±81 


8.28+0.12 


222030.68-004107.9 




17 


38±0.01 


8037±31 


8.28+0.07 


222551.65+001637.7 




18 


.95+0.02 


10640±94 


8.16+0.09 


231133.88+144119.9 




19 


.90+0.03 


39161±1369 


7.40±0.18 


233345.97-000843.0 




20 


.65+0.05 


11198±433 


8.35±0.29 



Helium fit. 
2 Blackb ody fit of SPSS pho tometry. 
3 lFarihietalJj2008l) 
4 IVanlandinghametalJ( 2005 ) 
5 IVennesetal.H2002l) 



combined white dwarf + blackbody model (WD+BB). This com- 
bined model represents a white dwarf with a cool (< 1000 K) disk. 
Starting at a temperature of 1000 K, the temperature of the sec- 
ondary blackbody was reduced in steps of 100 K until the combined 
WD+BB model matched the H and JC-band photometry. Here we 
assumed there is no excess in the JJ-band. This temperature can 
then be adopted as a upper limit as lower temperature disk combi- 
nations can also produce models which match the i\"-band photom- 
etry alone. NIR i\"-band spectroscopy or mid-infrared photometry 
would be needed to more accurately constrain the temperature of 
any debris disk containing system. 



5 FOLLOW-UP NIR PHOTOMETRIC OBSERVATIONS 

For a small number of sources, we obtained follow-up NIR 
photometry to independently measure their JHK magnitudes. 
SDSS J22255 1.65+001637.7 was observed on 2007 October 2 (at 
K s ) and 2007 October 3 with the 3.5m New Technology Tele- 
scope (NTT) and the Son-of-ISAAC (SOFI) Imager. Individual ex- 
posure times were 15 sees, 12 sees and 8 sees for the J, H and 
A's-bands respectively, with a total exposure time in each band of 
150 sees, 288 sees and 240 sees. The weather was clear with seeing 
of J — 1.3" on the first night and K s = 1.4" on the second night. 

SDSSJ003902.47-003000.3 was observed also with the 
NTT+SOFI on 2007 October 25 in J,H & K s . Individual ex- 
posure times were 15 sees, 12 sees and 12 sees for the J, H and 
K B -bands respectively, with a total exposure time in each band of 
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Spectral Type 

Figure 5. Absolute i^-band photometry verses spectral type for L and T 
dwarfs (Crosses: Reproduced from Vrba et al. 2004). Also plotted are the 
predicted absolute magnitudes of the candidate companions in this survey 
(Squares). 

75 sees, 144 sees and 144 sees. The weather was clear and photo- 
metric with seeing at K s = 0.7". 

WD 1318+005 was observed at the 3.9m Anglo-Australian 
Telescope (AAT) with the IRIS2 NIR imager on 2007 August 27. 
Individual exposures times were 5 sees in the A"-band with a total 
exposure time of 75 sees. The average seeing over the course of the 
observations was = 2.4". 

Standard nodding was used in all observations to ease te re- 
moval of NIR background. We reduced all the data automatically 
in the standard manner with the ORAC-DR package. The combined 
frames were photometrically calibrated using UKIDSS sources in 
each field. Aperture photometry was performed interactively using 
the Starlink Gaia package. 



6 CANDIDATE UKIDSS PHOTOMETRY 

For each candidate system with a putative stellar/substellar com- 
panion, we have estimated an upper limit on the projected orbital 
separation. This was achieved by measuring the FWHM of the sys- 
tem in the bandpass where the contributions of flux from each ob- 
ject would be roughly equal (usually the z' or J-band depending on 
the availability of the data), and equating this to an orbital sepera- 
tion at the estimated distance to the white dwarf. In the few cases 
where the candidate system has been resolved, the separation be- 
tween components has been measured directly and equated to a 
projected orbital seperation. 

For the few resolved systems, we have also calcu- 
lated the probability of a chance alignment (P a i) by refer- 
ring to the spatial number densities of M, L and T-dwarfs 
(Caballero, Burgasser & Klement 2008). The expected number of 
stars with the spectral type assigned to the putative companion is 
then calculated within a region of space defined by the errors on 
the white dwarf's distance. This then has to be multiplied by the to- 



tal sample size in order to take the complete volume of the survey 
into account. 

We discuss a number of examples in detail in the following 
sections. A complete summary of results is given in Table [7] and 
Figures inluded in Appendix lAl 

6.1 Examples I: White Dwarfs with Putative Low Mass 
Companions 

6. 1. 1 SDSSJ003902.47 -003000.3: A White Dwarf with an 
Unresolved Low Mass Stellar Companion 

SDSS J0039— 003's optical spectrum shows no obvious evidence 
for an M dwarf companion. EIS06 give the white dwarf physical 
parameters as T cB = 12, 314 ± 268K and log g = 7.39 ± 0.08 
from an automated fit to the optical spectrum. When compared to 
an atmospheric model at this temperature and gravity, a small ex- 
cess emission is visible in the SDSS z'-band (Figure lAH . with a 
clear strong excess emission in the UKIDSS YJHK-band pho- 
tometry. The predicted absolute K-band photometry of the poten- 
tial companion suggests a spectral type of M7. However, we found 
that a composite WD+L0 spectrum is the best fit to the UKIDSS 
photometry. 

We first identified this object in UKIDSS DR2 and so follow- 
up JH A'-band photometry was obtained using SOFI on the NTT in 
October 2007 (Section[5]l- The SOFI photometry confirms the pres- 
ence of a NIR excess and that the companion is likely of spectral 
type L0. The UKIDSS mergedClass statistic for SDSS J0039-003 
indicates the system is a point source, with a measured FWHM of 
0.73" in the J-band. This equates to a projected orbital separation 
of < 284 AU for the secondary at the estimated distance to the 
white dwarf (Table [8j, However, the low mass of the white dwarf 
(0.34±0.02A/q) suggests an accelerated evolution caused by a sec- 
ondary component, and so it is likely that the companion is close. 
If this is the case, a discrepancy in estimated spectral type of the 
secondary may arise from variability in the NIR photometry, due 
to night and day side heating of the companions atmosphere. NIR 
time series photometry is needed to confirm this scenario. 

6.1.2 SDSSJ012032.27-001351.1: A White Dwarf with a 
Resolved Low Mass Stellar Companion 

SDSSJ0120-001 was classified as a DA white dwarf by EIS06 
withTrff = 10,487±194Kandlogg = 7.91±0.23. However, the 
SDSS finder image shows a faint red background object which may 
have contaminated the optical spectrum to some extent (Figure |A3l 
An excess is clearly seen in the UKIDSS photometry, most likely 
due to this background source. Indeed, a closer inspection of the 
SDSS and UKIDSS images (Figure [6] & [7) reveals that the white 
dwarf has not been detected in the the NIR and hence it is only the 
flux from the nearby red object that is being measured. The two 
components are separated by ~ 2" equating to a projected orbital 
separation of ~ 910 AU. If the two are associated then the UKIDSS 
photometry is best matched by an M7 dwarf scaled to the estimated 
distance of the white dwarf. 

The expected numb er density of stars of spectral type M6- 
M8 V is 6.8 x 10~ 3 pc" 3 dCaballero et ai]2 008). Therefore, assum- 
ing an error of ±1 on the spectral type of the red object, the ex- 
pected number of stars of spectral type M6-M8 V in a volume of 
space defined by the errors of the white dwarf's distance (140 pc) 
and the projected separation is 5.8xl0~ J . Thus, this could be the 
widest white dwarf + M dwarf binary to be identified in this work. 




Figure 6. SDSS finding chart for SDSS J012032.27-001351.1. A sec- 
ondary object can be seen just to the south of the central white dwarf. The 
colours have been inverted from the original finding chart. 
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Figure 7. UKIDSS YJHK images of SDSS J012032.27-001351.1. Both 
objects from the SDSS image can be seen in the UKIDSS y-band, with 
only the southern most object appearing in the JHK images. 

Follow-up proper motion measurements are required to assess the 
association of the two objects, as a background giant star cannot be 
discounted at this stage. 

6.1.3 SDSS J013532. 98+144555. 8 = NLTT5306: A White 
Dwarf with an Unresolved Brown Dwarf Companion 

NLTT 5306 is classified by EIS06 as a DA white dwarf with T cS = 
8100±20 K and logg = 8.08±0.04. The SDSS optical photometry 
and spectrum (Figure lA4l show no sign of an excess due to a stellar 
companion. However, an excess is clearly visible in the UKIDSS 
photometry, with the predicted absolute it-band magnitude of the 
companion suggesting a spectral type of L6. A best fit is achieved 
with a WD+dL5 spectrum. At the white d warf age of 2.1 Gyr 
(Section[7}, the DU STY atmospheric models JChabrier et al.l2000l : 
Baraffe et al]|2002h predict a mass of 58±2Af Jllp . Therefore, it is 
likely SDSS 0135+144 harbours a brown dwarf companion. The 
UKIDSS mergedClass statistic for NLTT 5306 indicates the system 
is a point source, with a measured FWHM of 0.96" in the J-band. 
This equates to a projected orbital separation of < 57 AU for the 
secondary at the estimated distance to the white dwarf (Table[8}. 

NLTT 5 306 is known to have a high proper motion (^t a bs = 
0.194" vr' MLepine & Sharal2005n . The star was identified in both 
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Figure 8. POSS I and POSS II images of NLTT 5306, with the current and 
prior positions of the star respectively marked. Also shown is the current 
epoch finder image from SDSS, and a composite JHK image from the 
UKIDSS DR8. 

the first and second Palomar Sky Surveys (POSS I and POSS II), 
and a composite image was produced from the data (Figure [8j. 
No background object can be seen to be present along the path of 
NLTT 5306, increasing the likelihood that this star has a secondary 
low mass companion. 

6.1.4 SDSS J023247. 50-003909. 3: Contamination from a 
Foreground/Background Source 

SDSS 10232-003 is a DA white dwarf with EIS06 parameters of 
T eB = 9898 ± 161 K and logg = 8.04 ± 0.24. The white dwarf 
has i' — 20.11 ± 0.03 and so would be too faint to be detected 
by UKIDSS in the NIR. However, a secondary object was detected 
with a separation of 2.0" from SDSS J0232-003 (Figure s|9l&[T0l 
equating to a projected separation of 700 AU. At the estimated dis- 
tance to the white dwarf of 334 ± 83 pc, the if-band magnitude of 
the potential secondary indicates a spectral type of M8-L2. How- 
ever, it was found that a star of spectral type L3 was the best fit to 
the UKIDSS photometry. This difference in estimates could be due 
to the large uncertainties in the UKIDSS magnitudes and therefore 
distance estimate to the white dwarf. However, it is highly unlikely 
that an object of this spectral type would have been seen in the op- 
tical images. Therefore, it is more likely that this is a foreground 
or background object and so will be excluded from our statistical 
analysis. 

6.1.5 SDSSJ154431.47+060104.3: A White Dwarf with an 
Unresolved T Dwarf Companion 

SDSS 1544+060 shows no optical excess when compared to a 
blackbody model at T c g = 7946 ± 66 K, the temperature assigned 
by the EIS06 automated fit (Figure lA24t . An excess is clearly vis- 
ible in each of the UKIDSS bands with the predicted absolute K- 
band magnitude of the secondary suggesting a companion type of 
early T. Each data point matches a composite WD+T3 model par- 
ticularly well with the J-band photometry possibly indicating the 
detection of the CH4 absorption feature. If confirmed this would 
potentially be the first detection of a close WD+dT dwarf binary, 
and only the second detection of a WD+dT system. The UKIDSS 
mergedClass statistic for SDSSI1329+123 indicates the system is 
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Figure 9. SDSS finder image of SDSS J023247.50-003909.3. A secondary 
object can be seen just to the south west of the central white dwarf. The 
colours have been inverted from the original finding chart. 



Figure 11. /{"-band image of SDSS J22255 1.65+001637.7 taken with SOFI 
on the NTT. The star is in the centre and is clearly elongated, possibly indi- 
cating the system is becoming resolved. 




bital separation of < 435 AU for the secondary at the estimated 
distance to the white dwarf (Table [8}. 

Since SDSS J2225+001 is a good candidate for a rare, previ- 
ously unknown unresolved substellar companion to a white dwarf, 
independent JHK photometry was obtained at the 3.5m NTT at 
La Silla in October 2007 (Section[5j. Indeed, the SOFI JHK pho- 
tometry confirms the NIR excess and is a good fit to the WD+L7 
composite model. The SOFI A"-band image of SDSS J2225+001 
is shown in Figure [TT] and shows that the system may indeed be- 
come resolved with a higher resolution instrument. The FWHM of 
this image is 1.8" which equates to a projected orbital separation of 
< 350AU. 




Figure 10. UKIDSS JHK images of SDSS J023247.50-003909.3. The 
white dwarf has almost completely faded from view in the J-band image, 
with the red object completely dominating in the H and ii"-bands. 



a point source, with a measured FWHM of 2.12" in the J-band. 
This equates to a projected orbital separation of < 190 AU for the 
secondary at the estimated distance to the white dwarf (Table[8]l. If 
the system were in a close orbit the secondary would be the low- 
est object to have survived common envelope evolution. We now 
urgently require spectroscopic confirmation for this potential new 
type of WD binary system, with follow up radial velocity measure- 
ments to assess the separation of components. 



6.7.6 SDSS J222551. 65+001637.7: A White Dwarf with a 
Partially Resolved Brown Dwarf Companion 

SDSS J2225+001 shows no sign of an M dwarf in it's optical spec- 
trum and an atmospheric model with parameters T c s = 10640 ± 
94 K and \ogg = 8.16 ± 0.09 matches the SDSS photometry 
(Figure |A29t . A clear excess is then seen in the H and AT-band 
photometry with the predicted absolute A"-band magnitude of the 
secondary suggesting a spectral type of approximately L6. A com- 
posite WD+L7 model provides the best match in this case. The 
UKIDSS mergedClass statistic for SDSS J2225-001 suggest the 
system is partially resolved. However, there is not a clear separation 
of the components in the UKIDSS images. The measured FWHM 
of the star is 2.28" in the J-band which equates to a projected or- 



6.2 Examples II: Magnetic White Dwarfs 

We detected 19 hydrogen atmosphere magnetic white dwarfs 
(DAHs) in the cross correlation between the EIS06 and MS99, and 
the UKIDSS DR8 database. All of these were found to fit their 
predicted NIR photometry to within 3a of their UKIDSS observed 
values, with the exception of SDSS J121209.31+013627.7 (See be- 
low). In an attempt to increase this sample size, we performed a 
further cross correlation between the UKIDS S DR8 database an d 
the list of known magnetic white dwarfs from lKawka et alj J2007h . 
However, after accounting for duplicate entries between the cata- 
logues, this only increased the sample by a single DAH star. This 
star turned out to exhibit a NIR excess and we discuss it below. 



6.2.7 SDSS J '121209. 3 '7 +07 3627. 7 

SDSS J1212+013 was first identified as a binary through an optical 
spectrum and a single photometric J-band detection, limiting th e 
spectral type of the secondary to L5 or later ^Schmidt et al. 2005). 
It was initially thought to be a detached system until Debes et all 
(2006) detected strong cyclotron emission, an in dication of mass 
transfe r, in if-band time series p hotometry. Both lKoen & Maxted 
d2006l) and lBurleighetalJ ( 2006) detected optical variability (Pe- 
riod « 90 mins) consistent with that of a polar in a l ow accretion 
state. The NIR excess was successfully modelled by Fari hi et all 
( 2008), who modelled the system as a DAH + a brown dwarf com- 
panion of spectral type L8. Although the evidence suggests that 
SDSS J1212+013 is a CV in a low state, the NIR observations do 
not rule out a detached system with accretion occurring via a wind. 
This system shall therefore be included in our statistical analysis. 
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6.2.2 SDSS J 125044.42+ 154957. 4 = WD 1248+161 

WD 1248+161 is a PA H white dwarf identified by 
IVanlandinpham et al.l d2005h in the SDSS DR3. It has an es- 
timated effective temperature of T e g = 10, 000 K and a polar 
magnetic field strength B p =20 MG. The SDSS optical spectrum 
and photometry are a reasonable fit when compared to a blackbody 
model generated at the spectroscopic temperature of the white 
dwarf, with a clear excess in all four of the UKIDSS bands 
(Figure |Ai8l We found that all of the UKIDSS data could not be 
matched effectively with any composite model, with a discrepancy 
existing between the JH and the A-band photometry. This could 
be explained by the existence of strong cyclotron emission due 
to the magnetic field of the white dwarf, increasing the flux in 
the A-band. Therefore, the predicted absolute magnitude of the 
secondary was calculated using the /J-band, which suggests a 
spectral type of M8. A composite WD+M8 model matches the 
UKIDSS JH photometry. However, the J-band data must be 
interpreted with caution as it was flagged as 'noise' in the UKIDSS 
archive. We cannot rule out another source of this excess emission, 
such as the presence of a debris disk, and so NIR spectroscopy is 
required to discriminate between these possibilities. 

The UKIDSS mergedClass statistic for WDJ1248+161 in- 
dicates the system is a point source, with a measured FWHM of 
1.57" in the //-band. This equates to a projected orbital separation 
of < 245 AU for the putative secondary at the estimated distance 
to the white dwarf (Table [8j. 

6.3 Examples III: White Dwarfs with Debris Disks 

6.3.1 SDSS J 122859.92+ 104033.0 

SDSSJ1228+104 was first identified by iGiinsicke et alj J2006I) in 
DR6 of the SDSS, where they identified distinct emission lines 
of the Calcium 850-866 nm triplet. The lines profiles of the Ca 
emission triplet show a double-peaked morphology, indicative of 
a gaseous, rotating disk. The velocity of the Ca II lines peaks in- 
dicated an outer radius of t he disk of ~ 1.2_Rp. Th e star was later 
investigated in the NIR bv lBrinkworthetal.1 {2009), where it was 
discovered to have a large near- and mid-infrared excess, an indi- 
cation of the presence of a dusty component in addition to that of 
the gas. This was the first evidence for the coexistence of a both a 
gaseous and dusty debris disk around a white dwarf. 

We have recovered SDSS 11228+104 in this work as having a 
large A-band excess, as was expected from this system. 

6.3.2 SDSS J 132044.68+00 185 5.0 = WD 1318+005 

The SDSS optical spectrum of WD 1318+005 shows no evidence 
for a companion, and the UKIDSS JH photometry can be matched 
with a white dwarf model for T e « = 19, 649K and log g = 8.36 
(Figure lA20t . However, the UKIDSS A'-band measurement is 3. lo- 
in excess of this model, and can be matched by the addition of a 
600 K blackbody to the white dwarf model. The UKIDSS merged- 
Class statistic for WD 1318+005 indicates the system is a point 
source, which rules out contamination due to a background galaxy. 

We obtained NIR photometry of WD 1318+005 with IRIS2 
at the 3.9m Anglo-Australian Telescope (Section [5}- These data 
are consistent with both the A-band excess emission and the white 
dwarf model within errors (Figure |A20t , WD 1318+005 possibly 
possesses a warm debris disk, and if so then it is highly likely to 
be a DAZ white dwarf. Although it cannot be classified as a DAZ 



from the existing low resolution SDSS optical s pectrum, this star 
lies in the temperature regime of these objects. Zuckerma n et al.l 
(2003) suggest that « 25% of all DAs in this temperature range are 
metal-polluted DAZs, and higher resolution and higher signal-to- 
noise data are now urgently required to search for Call absorption. 

6.3.3 SDSS J 155720.77+091624.7 

SDSS J1557+091 shows no optical excess when fitted with an at- 
mospheric model with physical parameters of T e g — 21990 ± 
403 K and log g = 7.67 ± 0.06 (EIS06). An excess is clearly 
evident in the A-band photometry which is best matched by the 
addition of a 700 K blackbody, indicating a potential cool disk. 
However, we also note that there is a possible excess in the Y and 
J-band photometry which can be best matched by the addition of 
a companion with spectral type L4. This is within 2a of the H- 
band flux and so a sub-stellar companion, or indeed another source, 
cannot be discounted as a possibility at this stage. The UKIDSS 
mergedClass statistic for SDSS J1557+091 indicates the system is 
a probable point source (-2: 70% probability that the source is stel- 
lar), with a measured FWHM of 0.51" in the J-band. This equates 
to a projected orbital separation of < 250 AU for a secondary star 
at the estimated distance to the white dwarf (Table[8j. 



7 LOW MASS STAR OR BROWN DWARF? 

In order to assess whether or not each putative companion is a 
low mass main-sequence star or a brown dwarf an estimate of the 
mass of the secondary has been calculated. Firstly an age of the 
white dwarf was calculated by the addition of the white dwarf cool- 
ing age (Section [4. 2\ and the main-sequence lifetime of the pro- 
genitor sta r. This was es t imate d using the initial-final mass rela- 
tionship of Dobbie et al. (2006), which is valid for initial masses 
> 1.6Mn dKalirai et al.l 120081) . An approximate m ain-sequence 
lifetim e can then be calculated from the models of Girar di et al.l 
(2000). It should be noted that for white dwarfs where A/wd < 
0.5A/q it is highly likely that the star has evolved through mass 
transfer and for these stars an age can not be calculated through 
this method. However, this can be seen as further evidence for the 
existence of a secondary star. For these stars, a lower limit on the 
mass of the secondary is estimated by using the cooling age of the 
white dwarf. 

We then estimate a mass for the sec ondary by interpolat- 
ing the Lyon grou p atmospheric models dChabrier et al.1 l2000l : 
iBaraffe et al]|2002l) . given the age of the white dwarf and an esti- 
mate of effective temperature of the companion. We estimated these 
temperatures by comparison with observed M, L and T dwarfs 
dVrba et al.l l2004) and assuming an error of ±1 spectral type. For 
the WD+dM bi naries, we have estimated t he masses by referring 
to the models of IBaraffe & ChabrieJ H996). The results are given 
in Tableland plotted in Figure [121 

If we assume a value of 80Af Jup as a conservative cut-off mass 
for the stellar-substellar boundary, then 9-11 new DA white dwarf 
+ brown d warf binaries (1 DA +L8 has been spectroscopically con- 
firmed, see Steele et al. 2009). All of these have spectral types of 
M9 or later, with one candidate WD + T dwarf. 



8 A WIDE OR CLOSE BINARY? 

iFarihi et ail d2005l) suggest a bimodal distribution of separation for 
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Table 7. Summary of proposed sources for the near-infrared excesses identified in UKIDSS. For the resolved systems, the probability of a chance alignment 
is given. Further information on systems with a putative companion is given in Table 8 



SDSSJ WD Proposed Source of Excess Figure 

003902.47-003000.3 ■ Unresolved Companion IaTT 

003923.04+003534.7 - Unresolved Companion [A2l 

012032.27-001351.1 - Resolved Companion (P al = 5.8 x 10" 2 ) [A3] 

013532.98+144555.8 0132+145 Unresolved Companion |A4] 

032317.00-002612.7 - Foreground/Background Contamination 

033209.25+010722.1 Foreground/Background Contamination 

034221.55+005345.6 - Foreground/Background Contamination 

075327.05+244705.5 - Debris Disk |A5] 

082705.04+284402.7 0824+288 WD+dC+d2 (Farihi. Hoard & Wachter 20101 

085956.47+082607.5 - Unresolved Companion |A6] 

090759.59+053649.7 - Unresolved Companion |A7] 

092648.84+102828.8 - Unresolved Companion lA8l 

093821.34+342035.6 - Unresolved Companion |A9] 

100259.88+093950.0 - Unresolved Companion IaToI 

101642.93+044317.7 - Unresolved Companion IaTTI 

103220.19+011227.0 - Foreground/Background Contamination 

103448.92+005201.4 1032+011 Unresolved Companion IaTH 

103736.75+013912.2 - Unresolved Companion IaHI 

110826.47+092721.5 - Unresolved Companion lAlil 

113416.09+055227.2 - Unresolved Companion IA151 

115808.44-012312.9 1155-011 Unresolved Companion IA161 

120503.99+070031.0 Foreground/Background Contamination 

121209.31+013627.7 - DAH+L8 (Section 6.2.1) 

122625.93+002700.2 1223+007 Unresolved Companion 1X171 

122859.92+104033.0 - Debris Disk (Section 6.3.1) 

124455.15+040220.6 - Foreground/Background Contamination 

125044.42+154957.4 1248+161 Unresolved Companion (Magnetic WD) lAlU 

131955.04+015259.5 1317+021 Unresolved Companion [aT9| 

132044.68+001855.0 1318+005 Debris Disk lA20l 

132925.21+123025.4 - Unresolved Companion IA2T1 

134154.39+005600.8 1339+011 Unresolved Companion lA22l 

141448.25+021257.7 1412+024 Debris Disk IA231 

154431.47+060104.3 - Unresolved Companion lA24l 

155720.77+091624.7 - Debris Disk IA251 

162514.88+302610.8 Unresolved Companion lA26l 

220841.63-000514.5 - Resolved Companion (P al = 3.7 x 10" 1 ) lA27l 

222030.68-004107.9 Partially Resolved Companion (Steele et al. 2009) 

222551.65+001637.7 Partially Resolved Companion lA29l 

231133.88+144119.9 Foreground/Background Contamination 

233345.97-000843.0 - Unresolved Companion lA30l 



white dwarf + low mass secondary binaries. The UKIDSS merged- 
Class statistic may provide us with evidence which supports this 
conclusion. The white dwarf PHL 5038 was identified using this 
statistic as partially resolved, and was later fully resolved int o 
a wide orbiting (« 55 AU) WD+L8 binary dSteele etalj|2009h . 
Therefore, it may be possible to separate these stars into long and 
short period systems using this value. 

In order to investigate this possibility, the resolving power of 
UKIDSS must be determined for the binary candidates in this sur- 
vey. The UKIDSS photometry is measured using a 2" aperture in 
all bands. The projected orbital separation for this radius at dis- 
tances from 10-100pc are plotted in Figure [T3] The same is also 
plotted for the average FWHM of the binary candidates (~ 1.2"). 
We see that for the majority of stars it will not be possible to distin- 
guish between long and short period period systems if the system 
is classified as a point source, unless the binary has a particularly 
wide orbit. The maximum and projected orbital separations for all 
unresolved and resolved candidates in this survey are given in Ta- 
bled] 



9 LIMITS ON UNRESOLVED ULTRACOOL 
COMPANIONS TO DA WHITE DWARFS 

Atmospheric modelling successfully fitted 521 DA white dwarfs to 
within 3ct of both their UKIDSS H and K-band photometry. We 
will use these stars to estimate the limits on the spectral types of 
unresolved companions to DA white dwarfs and hence determine 
the sensitivity of UKIDSS to both L and T dwarfs in white dwarf 
binaries. The method used to calculate these limits is detailed be- 
low. 

In order to place a limit on the spectral type of the coolest 
detectable unresolved companion to each DA white dwarf we re- 
quired the absolute magnitu des for the spectral types of low mass 
objects ranging from L0-T8. |Patten et all 12006) lists known M, L 
and T dwarfs with 2MASS JHKs photometry ranging from spec- 
tral types M5-T8 and also has measured parallaxes for many of 
these. We converted the 2MASS JHKs photometry for each ob- 
ject to the U KIDSS JHK ph otometric system using the colour 
corrections of lCarpentej d200lb . Taking only those with a measured 
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Figure 12. Predicted masses of companions based on the Lyon group models. The dashed lines indicate constant age in Gyr. 
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Figure 13. Distance verses projected orbital separation for a 2" (solid) and 
1.2"(dashed) radius. 



parallax, and thus an estimated distance, the A-band photometry 
was scaled to 10 pc for each spectral type to obtain an estimated 
absolute A'-band magnitude. For multiple stars of the same spec- 
tral type we used an average of the magnitudes. For the few spec- 
tral types that had no parallax measurement we performed a linear 



interpolation to predict an absolute magnitude. The A-band mag- 
nitude was used for the reasons discussed in Section l4~3l 

For each apparently single DA white dwarf we calculated an 
expected A'-band magnitude (Section [4.2t . The observed A'-band 
magnitude was not utilised, as any of the white dwarfs could be 
in excess within the observed errors. This would have resulted 
in an over-estimate of the limiting spectral type for the compan- 
ion to these stars. We then calculated a 3a detection limit for 
each white dwarf by using the observed errors listed in UKIDSS 
DR8, and adding these to the predicted magnitudes. The distance to 
each white dwarf was calculated using the Bergeron models (Sec- 
tion !4.2b using the effective temperatures and surface gravities pro- 
duced from the automated fit of EIS06. For the single DA white 
dwarfs only found in MS99, temperature estimates were obtained 
from the literature or an automated fit to the stars optical photom- 
etry. For stars without a measured surface gravity, we assumed a 
logg = 8.0. This applied to only 11 white dwarfs in the sample, 
which is not a large enough number to create a distance bias that 
could later affect the statistics. Subsequently, the grid of M, L and T 
dwarf absolute magnitudes were scaled to the estimated distance of 
each white dwarf and added to the predicted white dwarf A'-band 
magnitudes until a match was made with the 3a detection limit. 
The spectral types of all the companion limits were summed (Ta- 
ble^ and plotted as a histogram and a cumulative histogram in 
Figure [T4l 

It should be noted that for the isolated white dwarf sample, the 
distances might be under-estimated if they are, for example, unre- 
solved double degenerate binaries. T hus, for this fra ction (sa 5% of 
all WDs are unresolved WD+WDs. rHolbersll2009h . the sensitivity 
to ultracool companions may have been overestimated. It is also 
important to spectroscopically confirm the candidate secondary in 
each candidate binary system for the same reason i.e. in case the 
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Table 8. Table of results for SDSS white dwarfs with NIR excess indicative of a companion: most likely companion spectral types and mergedClass stats 
indicating the possibility that a system is a widely separated or a close detached binary. 



SDSS J 


WD 


Spec. Types 


MergedClass 


M ffD (M s ) 


d(pc) 


Age (Gyr) 


McOMP (M Jup ) 


a(AU) 


10164? 0^4-044^17 7 

1\J 1 Ut^. "J yVrtTJ 11.1 






DA4-HM3 


_2 


40 4- 01 


uuu in ou 


> 0.01 


3£n + 1 ^»0 


< 530 


110826.47+092721.5 






DA+dM5 




1 18 + 02 


90 + 5 


2 6+ 01 


113 ± 23 


< 63 


162514 884-302610 8 






DA+riM5 




n 42 + 02 


(J80 + Qf) 

UUU _1_ i?U 


■> 0004 


113 ± 23 


< 545 


085956.47+082607.5 






DA+dM6 




0.36 ± 0.02 


400 ± 20 


> 0.1 


108 ± 18 


< 650 


113416.09+055227.2 






DA+dM6 


j 


0.57 ± 0.04 


580 ± 50 


1 6+°-l 

-0.4 


108 ± 18 


< 510 


012032.27—001351.1 






DA+dM7 




56 + 1 2 


410 ± 70 


9 9 + 5.8 


93 + 8 


~ 910 


103736.75+013912.2 






DA+dM7 




40 + 02 


330 + 10 


3 7+ 11 

°-'-0.8 


93 + 8 


< 330 


115808.44—012312.9 


1155- 


-011 


DA+dM7 


! 


«J 0.60 1 


280 ± 10 


« 1.8 1 


93 + 8 


< 270 


134154.39+005600.8 


1339+011 


DA+dM7 


+ 1 


0.58 + 0.02 


370 ± 20 


i.5+°i 

— 0.2 


93 + 8 


< 480 


132925.21+123025.4 






DA+dM7 




0.36 + 0.01 


210 + 10 


> 0.2 


93 + 8 


< 235 


233345.97-000843.0 






DA+dM7 


.] 


0.83 + 0.18 


390 ± 100 


1 1+ 1 - 2 

—0.4 


93 + 8 


< 350 


092648.84+102828.8 






DA+dM8 


.] 


0.68 ± 0.17 


480 ± 110 


—0.5 


85 + 5 


< 800 


125044.42+154957.4 


1248- 


1-161 


DAH+dM8 


.] 


k, 0.60 1 


150 ± 10 


> 0.6 


85 + 5 


< 245 


003923.04+003534.7 






DA+dM9 


-1 


0.60 + 0.06 


230 ± 20 


1 8+ ' 9 


80 + 3 


< 263 


003902.47-003000.3 






DA+dLO 




0.34 + 0.02 


390 ± 30 


> 0.3 


> 80 


< 484 


100259.88+093950.0 






DA+dLO 




0.57 + 0.03 


540 ± 30 


1 6 +0 - 5 


78 + 7 


< 420 


131955.04+015259.5 


1317- 


-021 


DA+dLl 




w 0.60 1 


290 ± 20 


» 1.3 1 


71 + 5 


< 130 


220841.63 -000514.5 






DA+dLl 


+ 1 


0.78 + 0.08 


170 ± 20 


i c+0.6 
i '°-0.3 


71 + 5 


290 


090759.59+053649.7 






DA+dL4 




0.54 + 0.02 


360 ± 20 


9 1+0.5 
" s - L -0.4 


58 + 2 


< 600 


013532.98+144555.8 


0132- 


-145 


DA+dL5 




0.65 + 0.03 


60 ± 10 


o 1+0.3 
Z - 1 -0.2 


58 + 2 


< 57 


093821.34+342035.6 






DA+dL5 




0.61 + 0.08 


160 ± 20 


2 2+ 1 - 4 
z ' z -0.7 


59 + 5 


< 154 


103448.92+005201.4 


1032- 


-011 


DA+dL5 




0.68 + 0.10 


200 ± 30 


1 4+0-9 
± '*-0.4 


55 + 4 


< 150 


222551.65+001637.7 






DA+dL7 


+ 1 


0.70 + 0.06 


190 ± 20 


i 9+0.4 
^-0.3 


47 + 3 


< 350 


121209.31+013627.7 






DAH+dL8 




rj 0.60 2 


150 2 


1-5 2 


50 2 


0.003 


222030.68-004107.9 






DA+dL8 


+ 1 


0.77 + 0.04 


60 ± 10 


9 1+0.4 
z ' 1 -0.3 


49 + 3 


55 


154431.47+060104.3 






DA+dT3 




0.85 + 0.09 


90 ± 10 


2 6+ ' 8 
Z -°-0.8 


46 + 2 


< 190 



Surface gravity estima ted from ma s s distr ibution of white dwarfs. 
2 Farihi et alj J2008b 




Figure 14. Left: Distribution of limits on unresolved companions to the 521 "single" DA white dwarfs (with well determined parameters) detected in UKIDSS 
DR8. Right: Plotted as a cumulative distribution. 
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Table 9. Number of each spectral type obtained as limits to the 'single' DA 
white dwarf population in UKIDSS DR8. 



Spectral Type 


No. 


Spectral Type 


No. 


Spectral Type 


No. 






LO 


1 


TO 


56 






LI 


7 


Tl 


41 






L2 


4 


T2 


48 






L3 


8 


T3 


33 


<M5 





L4 


21 


T4 


33 


M5 


1 


L5 


33 


T5 


26 


M6 


1 


L6 


21 


T6 


16 


M7 


2 


L7 


37 


T7 


15 


M8 


2 


L8 


48 


T8 


8 


M9 


1 


L9 


47 


>T8 


11 



distance to the white dwarf has been under-estimated, and thus the 
temperature and mass of the companion are under-estimated. 



10 SURVEY SENSITIVITY TO UNRESOLVED 
COMPANIONS 

Figure [14] shows that the sensitivity of this survey for the isolated 
DA white dwarf sample peaks roughly at the L-T borderline, at ap- 
proximately L8-T2. Figure[l4]also demonstrates that the sensitivity 
function then falls off as the limits approach a spectral type of late 
T, with only 43% of the sample sensitive to the detection of a com- 
panion earlier than or equal to spectral type TO. In contrast 99% of 
the sample is sensitive to the detection of a companion of spectral 
type earlier than or equal to LO, dropping to 5 1 % by L9. Figure [TBI 
shows a histogram of the distances to these single white dwarfs. 
98% of the sample lies within 500 pc of the Sun. 

In considering the total survey sensitivity, we must take the the 
274 DA+dM binaries that fell onto our model grid into account. 
These must be included in the final sample for statistical analysis 
as many would have been sensitive to the detection of L and T-type 
companions, if the M dwarf companion were to be removed. How- 
ever, the presence of the M dwarf companions may bias this popu- 
lation towards greater distances, as the secondary is far brighter in 
the near-infrared. Therefore, the fraction of DA+dM binaries that 
lie within 500 pc must be determined so a direct comparison can 
be made against the isolated white dwarf population. The distances 
to these were estimated using the effective temperatures and sur- 
face gravities assigned by EIS06. We plot this data as a histogram 
in Figure[l5] with 93% of the DA+Ms lying within 500 pc. There- 
fore, there is not a significant difference between the isolated and 
white dwarf + M dwarf distance distributions. 



11 BINARY STATISTICS FOR WHITE DWARFS WITH 
UNRESOLVED SUBSTELLAR COMPANIONS 

In this section we are going to estimate the fractions of DA white 
dwarfs with unresolved L dwarf, T dwarf and brown dwarf com- 
panions using the results of this survey. 

We assume that the white dwarf + L and T dwarf candidates, 
and the candidate disk systems, follow the same distance distribu- 
tion as the single white dwarfs, as shown in the previous section. 
The new candidate white dwarf + M dwarf systems can be added 
to the previously identified WD+dM binary sample. Therefore, the 
complete sample size of DA white dwarfs at any distance, with or 
without a near-infrared excess is 794. The number of these stars 



that lie within 500 pc is 98% of the total sample size giving 778 
DA white dwarfs. 

We can now calculate an upper limit for the unresolved L 
dwarf, T dwarf and brown dwarf companion fractions to DA white 
dwarfs in UKIDSS. This will be done firstly for the DA + L dwarfs, 
followed by the DA + T dwarfs and then finally for the DA + 
brown dwarf population. 

The candidate L and T dwarfs have been identified through 
NIR photometry with only PHL5038 and SDSS J1212+013 con- 
firmed as harbouring brown dwa rf companions. PH L5038B has 
been confirmed spectroscopically JSteele et al.ll2009h . and the sub- 
stellar status of secondary to S DSS 1212+013 is supported by a 
number of convincing arguments jFarihi et al.l2008h . However, due 
to its high proper motion, we have determined that the source of the 
near-infrared excess to NLTT5306 is directly associated with the 
star, and not due to a background object (Section[6jT3}- Therefore, 
we will assume that at least 3 of the detections in this survey are 
confirmed white dwarf + brown dwarf binary systems containing 
L dwarf secondaries. These confirmed brown dwarfs can be used 
to set a lower limit on the unresolved L dwarf and brown dwarf 
companion fractions. 

Table 1 101 shows the percentage of the sample sensitive to the 
detection of L dwarfs equal to or earlier than the specified spec- 
tral type, and the number of each spectral type potentially detected 
in this survey. In order to calculate the white dwarf + L dwarf bi- 
nary fraction, we are going to determine the effective number of 
detections. This value represents the number of each spectral type 
detected if the survey were 100% efficient at detecting companions 
of all spectral types. The effective number of detections can then be 
calculated by dividing the actual number of detections by the sen- 
sitivity to each spectral type. For example, if the survey sensitivity 
to companions of spectral types earlier or equal to L0 was 50%, 
and 1 L0 was detected, then the effective number of detections is 2. 
This is then summed for all spectral types to give the total effective 
number of L dwarfs detected: 

2 2 1 3 

UcB ~ 099 + 098 + 097 + 096 + 095 + 091 + 084 



0.80 0.73 0.64 



The L dwarf companion fraction is then the effective number of de- 
tections divided by the sample total. An error can be estimated by 
taking the square root of the reciprocal number of effective detec- 
tions. Thus, an upper limit to the fraction of DA white dwarfs with 
unresolved L dwarf companions is /wD+dL ^ 1.6+0.5%. A lower 
limit can be estimated by assuming only 3 candidates (PHL 5038, 
SDSS 1212+013 and NLTT 5306) are real. This gives a final range 
Of 0.4 < /WD+dL < 2.1%. 

Table [10] also gives the percentage of the sample sensitive to 
the detection of T dwarfs equal to or earlier than the specified spec- 
tral type, and the number of each spectral type potentially detected 
in this survey. An upper limit to the fraction of white dwarf with 
T dwarf companions is then calculated as detailed previously giv- 
ing /wD+dT ^ 0.8%. Although only one candidate T dwarf has 
been detected, an average over the entire T dwarf spectral range 
is a good start under these circumstances but ultimately, a more 
complex Bayesian analysis would yield a better result. Additional 
candidates from future surveys with greater sensitivity and a larger 
search field will greatly improve the statistics. 
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Figure 15. Left: Cumulative distribution of distances to the isolated DA white dwarfs in this work. Right: Cumulative distribution of distances to the DA+dM 
binaries in this work. 



In order to determine the fraction of white dwarfs with brown 
dwarf companions, we must first estimate the spectral type where 
the secondary becomes truly substellar. This depends on the age of 
the secondary, which can be estimated using the total age (main- 
sequence lifetime + cooling age) of the white dwarf primary. The 
average age of the white dwarf sample is 1.9±0.7Gyr. After ~ 
0.1 Gyr brown dwarfs cool very slowly, so the average sample age 
can be used to estimate a spectral type where the stellar/substellar 
borderline occurs. Taking the upper limit on brown dwarf mass as 
80M Jup , then an upper limit for the brown dwarf spectral type can 
be taken as L0 (Figure [Lit. 

An upper limit to the fraction of white dwarf with brown 
dwarf companions is then calculated as detailed previously giv- 
ing /wd+bd ^ 2.1 ± 0.5%. A lower limit can be estimated 
by assuming only 3 candidates (PHL5038, SDSS 1212+013 and 
NLTT 5306) are real. This gives a final range of 0.4 ^ /wd+bd ^ 
2.6%. 

Although these statistics are not particularly robust, they are 
suitable as a first approximation using the small numbers available. 
In order to improve upon these numbers, the sample needs to be 
enlarged. In the first instance, this will be done by UKIDSS, which 
is set to be completed by 2012. Looking ahead, future infrared sur- 
veys with VISTA and WISE will add to the white dwarf + ultracool 
companion sample significantly. 

It should also be noted that the calculations presented here 
underestimates the total fraction of white dwarfs with ultracool 
companions, as this search does not include widely orbiting, re- 
solved com panions. Searches for these systems are already being 
undertaken (Da y-Jones et al. 2008). Due to the aperture size used in 
analysing UKIDSS photometry (2"), we have found a small num- 
ber of such systems as part of this work. This is due to the merged- 
Class statistic which, as noted previously, provides an indication of 
whether a source is a point or extended source (Section[8]l. In the 
first instance this value can be used to remove contaminating galax- 
ies from a search, but it can also be an indication that a system is 
partially resolved (e.g. PHL 5038). Thus, for this survey wide com- 



panions have only been discovered if they happen to be partially 
resolved in the UKIDSS database. 



12 CORRECTING FOR COMPLETENESS 

Some of the white dwarfs with a near-infrared excess may not have 
been detected in UKIDSS if the companion were to be removed. 
We can correct for this bias by removing those stars with predicted 
AT-band magnitudes which are outside the 5<r detection limits of the 
survey. We can then repeat the statistics using this new, magnitude 
limited sample. After applying this constraint, we are left with 413 
"single" DA white dwarfs and 151 previousley identified DA+dM 
binaries, giving a total sample size of 564. Only 3 of the white 
dwarf + brown dwarf candidates (SDSS J0135+144, PHL 5038 
and SDSS J 1544+060) are within the magnitude limited sample, 
giving the following binary fractions; /wD+dL ^ 0.4 ± 0.3%, 
/wD+dT > 0.2% and /wd+bd > 0.5 ± 0.3. 

13 COMPARISON TO PREVIOUS ESTIMATES 

13.1 The White Dwarf + L dwarf Binary Fraction 

iFarihi et all d2005h (hereafter FBZ05) estimated the white dwarf + 
L dwarf binary fraction as /wD+dL ^ 0.5% from a MR study of 
371 white dwarfs with 2MASS photometry. FBZ05 estimated their 
average sensitive to unresolved companions to be L8, the same as 
the peak of the sensitivity in the UKIDSS data. In order for the L 
dwarf fraction from UKIDSS (0.4 < /wD+dL < 2.1%) to agree 
with the FBZ05 result either (1) All of the remaining photometric 
candidates presented here are not WD+dL binaries i.e. the lower 
limit of the calculated L dwarf fraction is more likely to be correct, 
or (2) FBZ05 have under-estimated the fraction of white dwarfs 
with L dwarfs. As the first candidate for which spectroscopy was 
obtained (PHL 5038) was confirmed to be real, and in other cases 
the photometric classification remains quite convincing, then ex- 
planation (1) seems initially less likely than (2). To investigate this 
further, we will perform a detailed examination as to how FBZ05 
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Table 10. Actual number and effective number of detections for DA white dwarf + dL/dT binaries in the UKIDSS DR8 sample based on sensitivity estimates. 
No candidate companions later than spectral type T3 were detected. 



Spectral 


No. Isolated WDs 


% Detectable of 


Actual No. 


Effective No. 


Sample % 


Type 


Sensitive to Spectral 


Spectral Type 


Candidates 


Detected 






Type or Earlier 


or Earlier 


Detected 






LO 


514 


99% 


2 


2.0 


0.3% 


LI 


513 


98% 


2 


2.0 


0.3% 


L2 


506 


97% 





0.0 


0.0% 


L3 


502 


96% 





0.0 


0.0% 


L4 


494 


95% 


1 


1.1 


0.1% 


L5 


473 


91% 


3 


3.3 


0.4% 


L6 


440 


84% 





0.0 


0.0% 


L7 


419 


80% 


1 


1.3 


0.2% 


L8 


382 


73% 


2 


2.7 


0.3% 


L9 


334 


64% 





0.0 


0.0% 


TO 


287 


55% 





0.0 


0.0% 


Tl 


231 


44% 





0.0 


0.0% 


T2 


190 


37% 





0.0 


0.0% 


T3 


142 


27% 


1 


3.7 


0.5% 


T4 


109 


21% 





0.0 


0.0% 


T5 


76 


15% 





0.0 


0.0% 


T6 


50 


10% 





0.0 


0.0% 


T7 


34 


7% 





0.0 


0.0% 


T8 


19 


4% 





0.0 


0.0% 



arrived at their L dwarf companion fraction. We will then follow 
the same method for the UKIDSS sample, and make a direct com- 
parison to the detailed sensitivity analysis discussed previously. 

It should be noted that the FBZ05 sample includes white 
dwarfs not of spectral type DA. However, these are a minority in the 
sample, and their removal would not significantly affect the fraction 
calculated by FBZ05. 

Figures [16] and [TTJ show the distance distributions of both 
surveys, with the mean distances plotted (FBZ05 = 57 pc and 
UKIDSS = 161 pc). FBZ05 calculate the detection limit of their 
survey by determining the "average" white dwarf from their sample 
(found by calculating the average temperature, T e g — 13000 K, of 
their sample and assuming logg = 8.0), and placing an upper limit 
on any putative companion's spectral type assuming this "average" 
white dwarf is situated at the mean sample distance. FBZ05's "av- 
erage" white dwarf has absolute magnitudes Mh = 11.8 mags 
and Mk b = 11.9 mags. These magnitudes are beyond the detec- 
tion limits of the 2MASS survey (Mh ~ 11.6 mags and K s ~ 
10.8 mags), and so FBZ05 set the detection limit by assuming a 
companion has a magnitude equal to the difference between the 
2MASS survey limits and this "average" white dwarf. This gives a 
limit of L8 and L4 in the 2MASS H and AVbands respectively. 

Following the FBZ05 method, an "average" white dwarf must 
be analysed using the UKIDSS sample, for which the average 
temperature is T c g — 14000 K. Assuming a surface gravity of 
loggr = 8.0 this equates to 2MASS absolute magnitudes of Mh = 
11.7 mags and Mk s = 11.8 mags. Using the average distance of 
161 pc, the 5a limits of the UKIDSS survey are Mh = 12.8 mags 
and Mk s = 12.2 mags. 

It should be noted at this point that these are fainter than the 
corresponding magnitudes of the "average" white dwarf, differing 
from the FBZ05 analysis. In FBZ05 the "average" white dwarf is 
in fact brighter in magnitude than the 2MASS limits, allowing for a 
limit to be set on the minimum luminosity of a putative companion 
by using the difference between those two values. As the opposite 
case is true for UKIDSS, it would make sense to use the addition 
of a typical 3a error as the limiting magnitude. An interesting point 



can be made here; if FBZ05 had added the 3a error at the 2MASS 
limit to their initial estimate, then the "average" companion would 
have been brighter, and therefore of an earlier spectral type. 

At 161 pc the "average" UKIDSS white dwarf has magnitudes 
H — 17.7 mags and K = 17.8 mags, which have typical errors of 
±0.05 and ±0.1 respectively. Assuming a companion is detected 
at the > 3<t level, this corresponds to a companion with apparent 
magnitudes of H — 20.0 mags and K — 18.8 mags. This sets a 
limiting spectral type for the UKIDSS survey of later than T9 in 
both bands. 

Figure [14] shows that it is indeed the case that the UKIDSS 
survey is sometimes sensitive to late (>T8) T dwarf companions. 
However, the sensitivity analysis for each white dwarf individually 
shows that this accounts for < 2% of the entire sample. Clearly it 
is not the case that this survey is capable of detecting late T dwarf 
companions around every white dwarf. Indeed, we have shown in 
Section[l0]that the sensitivity of this survey peaks at the L-T border. 
The UKIDSS "average" white dwarf is clearly not representative of 
either the mean or the median star in the sample. Indeed, it is only 
representative of the extreme of the survey sensitivity and therefore 
not representative at all. Therefore, we suspect that FBZ05's "av- 
erage" white dwarf is unlikely to represent their mean or median 
star either. For example, if the distance distribution to each sample 
is inspected, 35% and 33% of white dwarfs lie beyond the average 
distance from FBZ05 and UKIDSS respectively. Hence, it is likely 
that at least that fraction are not in fact as sensitive to low mass 
companions as the "average" star. 

FBZ05 state "NIR excess detection requires photometric ac- 
curacy, not deep imaging". In fact, as UKIDSS is ~3 magnitudes 
deeper than 2MASS, the UKIDSS photometric errors are much 
smaller at the fainter end of the 2MASS survey where many of the 
white dwarfs lie. In 2MASS the photometric signal-to-noise ratio 
is > 10 for Ks < 14.3 mags, so a typical error is 0.1 mags. In 
the equivalent UKIDSS bandpass, the typical error is 0.005 mags. 
Clearly, UKIDSS will provide much better accuracy at the 2MASS 
limits, and continue this trend until the limits of the survey are 



18 P. R. Steele et al. 



1 20 




500 



Distance (pc) 

Figure 16. Distance distribution of the white dwarfs presented in Farihi et 
al. (2005). The mean distance is shown as a solid line. 
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Figure 17. Distance distribution of the single DA white dwarfs present in 
the UKIDSS sample. The mean distance is shown as a solid line. 



reached. Indeed, it may be the case that 2MASS is not as sensitive 
as the "average" white dwarf would suggest. 

In conclusion, the use of an "average" white dwarf to place 
sensitivity limits on a sample may be the problem. The detection 
of an ultracool companion requires that the companion exceed the 
white dwarf's luminosity by 3a at an infrared wavelength. The lu- 
minosity of a star is proportional to R and T 4 , hence an "average" 
white dwarf determined according to temperature is not represen- 
tative of the sample. Unfortunately, FBZ05 do not detail the pa- 
rameters needed for each white dwarf (Effective temperature and 
distance) to perform a similar analysis to that presented for the 



UKIDSS data. The reasons for the surprisingly similar sensitivities 
of both the 2MASS and UKIDSS surveys cannot be fully under- 
stood until such an analysis can be performed. 



13.2 The Main-sequence Star + Brown Dwarf Binary 
Fraction 

iMcCarthv & Zuckermatj J2004h presented a coronographic NIR 
search for substellar companions to nearby main-sequence stars. 
They found that the frequency of brown dwarf companions to 
G, K and M stars orbiting at a distance between 75-300 AU to 
be 1±1%. They also calculated a companion fraction for mas- 
sive brown dwarfs orbiting be tween 120-1200 AU to be 0.7±0.7%. 
iGrether & Lineweaver] d2006t) analysed the relative fraction of close 
(P< 5 yr) stellar, brown dwarf and planetary companions to nearby 
sun-like stars and found that close brown dwarf companions occur 
at a frequency of < 1%. Therefore, the brown dwarf companion 
fraction to white dwarfs (0.4 ^ /wd+bd ?S 2.6%) and hence, their 
progen itors, from the UKIDSS survey agr ees with the estimates 
of bothlMcCarthv & ZuckermarJ (2004) and Grether & Lineweaverl 
d2006h within errors. 

However, there exists a bias among white dwarfs in that they 
are generally descended from stars more massive than the Sun. It 
could easily be the case that such stars have a greater frequency of 
brown dwarf companions than solar-type stars. 

The brown dwarf desert is a noted deficit in the frequency of 
brown dwarf compani ons relative to the freq uency of less massive 
planetary companions ( Marcv & Butler 2000) or relative to the fre- 
quency of more mass i ve ste llar companions to Solar-like hosts. 
I Grether & Lineweaver] (20061) verified the existence of this phe- 
nomenon, with the desert extending to < 3AU. Determining the 
separations of the WD+BD binaries in the UKIDSS sample will 
allow additional investigation of the brown dwarf desert. It may 
even be possible to determine the companion fraction at different 
orbital radii. 



14 SUMMARY 

We have cross correla ted the UKIDSS LAS DR8 with both the 
lEisenstein et alj (2006) and the McCook & Sion online catalogues 
of spectroscopically identified white dwarfs, in order to search for 
stars displaying near-infrared excesses indicative of low mass com- 
panions or debris disks. The sample included 811 hydrogen at- 
mosphere DA white dwarfs with 1500 < T e g < 100, 000 K, 
6.5 <logg < 9.5, and both H and A"-band photometry present 
in the UKIDSS archive. 

We identified 313 white dwarfs with near-infrared excesses, 
with 274 of these recovering previously identified binaries. This 
left 39 candidate white dwarfs for further analysis. Twenty seven 
of these had multiple excesses indicative of a low mass compan- 
ion, with 9-11 of these having a predicted mass in the range as- 
sociated with brown dwarfs, and a further 14-16 likely very low 
mass stellar companions. Seven stars in the sample showed ev- 
idence of contamination by a foreground or background object. 
The remaining 3 each had a A'-band excess indicative of a debris 
disk including the already discovered DAd SDSS J1228+104. Two 
magnetic white dwarfs were found with a near-infrared excess; (1) 
SDSSJ1212+01 - a previously identified DAH+dL8 binary, and 
(2) SDSS J1250+ 154 - potentially explained by the presence of an 
M8 companion and additional cyclotron emission. 
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We have placed limits on the latest spectral type of unre- 
solved companion detectable around each of the apparently sin- 
gle DA white dwarfs detected in the UKIDSS database. We found 
that this survey has a sensitivity of » 99% to companions of 
spectral type LO or earlier and « 43% to companions of spec- 
tral type TO or earlier. These sensitivities were used in conjunc- 
tion with the photometrically detected L and T dwarf companions 
to estimate the unresolved brown dwarf companion fraction to DA 
white dwarfs in UKIDSS as 0.4 ^ /wd+bd 2.6%. Correcting 
for the completeness of the survey gives a final binary fraction of 
/wd+bd > 0.5 ± 0.3%. 

The W D+dL binary fractio n in UKIDSS is significantly higher 
than that of iFarihi et al.l d2005t). One po ssible explanation for this 
discrepancy is that IFarihi et al.l d2005f) have used an "average" 
white dwarf to set the limits of their 2MASS survey. By utilising 
the same method for the UKIDSS white dwarfs, we have shown that 
this "average" white dwarf may not be representative of the mean 
survey sensitivity. A more thorough analysis of the 2MASS survey 
sensitivity is required, using the same method employed here, be- 
fore the differences can be fully understood. This will require each 
white dwarf to be analysed individually in order to place a limit on 
the coolest detectable companion around each star. 

We found that the fraction of ultracool companions to white 
dwarfs was found to be consistent with estimates of the main- 
sequence + brown dwarf binary fraction. 



15 NOTES 

An updated verion of Pierre Bergeron's synthetic 
model grid was made available in February 2011 
dTremblavTB ergeron & Gianninas! 1201 ll) . These revisions should 
have little effect on the magnitudes and distances calculate in this 
work, so will be used for the analysis of future UKIDSS data 
releases. 

Before submission the authors became aware of the study of 
Girven et al. (2011, submitted) who have carried out an indepen- 
dent study to identify photometric white dwarf candidates in SDSS 
DR7, and also search for infrared excesses. A comparison of results 
shows each study to be broadly consistent. 



16 ACKNOWLEDGEMENTS 

PRS is supported by RoPACs, a Marie Curie Initial Training Net- 
work funded by the European Commissions's Seventh Frame- 
work Programme. MRB acknowledges receipt of an STFC Ad- 
vanced Fellowship. Balmer lines in the models were calculated 
with the modified Stark broadening profiles of Tremblay & Berg- 
eron (2009), kindly made available by the authors. 



REFERENCES 

Baraffe I., Chabrier G., 1996, ApJL, 461, L51+ 
Baraffe I., Chabrier G., Allard F, Hauschildt P. H., 2002, A&A, 
382, 563 

Becklin E. E., Farihi J., Jura M., Song I., Weinberger A. J., Zuck- 

erman B., 2005, ApJL, 632, L119 
Becklin E. E., Zuckerman B., 1988, Nature, 336, 656 
Bergeron P., Saffer R. A., Liebert J., 1992, ApJ, 394, 228 
Bergeron P., Wesemael F, Beauchamp A., 1995, PASP, 107, 1047 



Brinkworth C. S., Gansicke B. T, Marsh T. R., Hoard D. W., Tap- 
pert C, 2009, ApJ, 696, 1402 

Burleigh M. R., et al., 2006, MNRAS, 373, 1416 

Burleigh M. R., Hogan E., Dobbie P. D., Napiwotzki R., Maxted 
P. F. L., 2006, MNRAS, 373, L55 

Caballero J. A., Burgasser A. J., Klement R., 2008, A&A, 488, 
181 

Carpenter J. M., 2001, AJ, 121, 2851 
Casali M., et al, 2007, A&A, 467, 777 

Chabrier G., Baraffe I., Allard E, Hauschildt P., 2000, ApJ, 542, 
464 

Chiu K, Fan X., Leggett S. K, Golimowski D. A., Zheng W, 
Geballe T. R., Schneider D. P., Brinkmann J., 2006, AJ, 131, 
2722 

Cushing M. C, Rayner J. T, Vacca W. D., 2005, ApJ, 623, 1115 

Day-Jones A., Pinfield D. J., Jones H. R. A., Napiwotzki R., Burn- 
ingham B., Jenkins J. S., UKIDSS Cool Dwarf Science Working 
Group 2008, in American Astronomical Society Meeting Ab- 
stracts Vol. 21 1 of American Astronomical Society Meeting Ab- 
stracts, Hunting For Wild Brown Dwarf Companions To White 
Dwarfs In UKIDSS And SDSS. pp 103.30-+ 

Day-Jones A. C, Pinfield D. J., Ruiz M. T, Beaumont H., Burn- 
ingham B., Gallardo J., Gianninas A., Bergeron P., Napiwotzki 
R., Jenkins J. S., Zhang Z. H., Murray D. N., Catalan S., Gomes 
J., 2011, MNRAS, 410, 705 

Debes J. H., Lopez-Morales M., Bonanos A. Z., Weinberger A. J., 
2006, ApJL, 647, L147 

Dobbie P. D., Burleigh M. R., Levan A. J., Barstow M. A., Napi- 
wotzki R., Holberg J. B., Hubeny I., Howell S. B., 2005, MN- 
RAS, 357, 1049 

Dobbie P. D., Napiwotzki R., Burleigh M. R., Barstow M. A., 

Boyce D. D., Casewell S. L., Jameson R. F, Hubeny I., Fontaine 

G„ 2006, MNRAS, 369, 383 
Dye S., et al„ 2006, MNRAS, 372, 1227 
Eisenstein D. J., et al., 2006, ApJS, 167, 40 
Farihi J., Becklin E. E., Zuckerman B., 2005, ApJS, 161, 394 
Farihi J., Burleigh M. R., Hoard D. W., 2008, ApJ, 674, 421 
Farihi J„ Christopher M„ 2004, AJ, 128, 1868 
Farihi J., Hoard D. W., Wachter S., 2006, ApJ, 646, 480 
Farihi J., Hoard D. W., Wachter S., 2010, ApJS, 190, 275 
Gansicke B. T, Marsh T. R., Southworth J., Rebassa-Mansergas 

A., 2006, Science, 314, 1908 
Girardi L., Bressan A., Bertelli G., Chiosi C, 2000, AAS, 141, 

371 

Golimowski D. A., et al., 2004, AJ, 128, 1733 

Green P. J., Ali B., Napiwotzki R., 2000, ApJ, 540, 992 

Grether D., Lineweaver C. H., 2006, ApJ, 640, 1051 

Hewett P. C, Warren S. J., Leggett S. K, Hodgkin S. T, 2006, 

MNRAS, 367, 454 
Hoard D. W., Wachter S., Sturch L. K., Widhalm A. M., Weiler 

K. P., Pretorius M. L., Wellhouse J. W., Gibiansky M., 2007, AJ, 

134, 26 

Holberg J. B., 2009, Journal of Physics Conference Series, 172, 
012022 

Holberg J. B., Bergeron P., 2006, AJ, 132, 1221 

Holberg J. B., Magargal K, 2005, in 14th European Workshop 

on White Dwarfs, ed. D. Koester & S. Moehler, ASP Conf. Ser. 

(San Francisco ASP) Vol. 334, . p. 419 
Hubeny I., 1988, Computer Physics Communications, 52, 103 
Hubeny I., Lanz T, 1995, ApJ, 439, 875 
JuraM., 2003, ApJ, 584, L91 

Jura M., Farihi J., Zuckerman B., 2007, ApJ, 663, 1285 



20 P. R. Steele et al. 



Kalirai J. S., Hansen B. M. S., Kelson D. D., Reitzel D. B., Rich 

R. M, Richer H. B., 2008, ApJ, 676, 594 
Kawka A., Vennes S., Schmidt G. D., Wickramasinghe D. T., 

Koch R., 2007, ApJ, 654, 499 
Kilic M., Redfield S., 2007, ApJ, 660, 641 

Kilic M, von Hippel T., Leggett S. K., Winget D. E., 2005, ApJ, 
632, LI 15 

Kilic M., von Hippel T., Leggett S. K., Winget D. E., 2006, ApJ, 
646, 474 

Knapp G. R., et al., 2004, AJ, 127, 3553 

Knutson H. A., et al., 2007, Nature, 447, 183 

Koen C, Maxted P. F. L., 2006, MNRAS, 371, 1675 

Koester D., 2008, ArXiv e-prints 

Lawrence A., et al., 2007, MNRAS, 379, 1599 

Lepine S., Shara M. M., 2005, AJ, 129, 1483 

Liebert J., et al., 2005, AJ, 129, 2376 

Marcy G. W., Butler R. P., 2000, PASP, 112, 137 

Maxted P. F. L., Napiwotzki R., Dobbie P. D., Burleigh M. R., 

2006, Nature, 442, 543 

McCarthy C., Zuckerman B., 2004, AJ, 127, 2871 

McCook G. P., Sion E. M., 1999, ApJS, 121, 1 

Patten B. M., Stauffer J. R., Burrows A., Marengo M., Hora J. L., 
Luhman K. L., Sonnett S. M., Henry T. J., Raghavan D., Megeath 
S. T., Liebert J., Fazio G. G, 2006, ApJ, 651, 502 

Probst R. G„ 1983, ApJS, 53, 335 

Reach W. T., Kuchner M. J., von Hippel T., Burrows A., Kilic M., 
Mullally R, Winget D. E., 2005, ApJ, 635, L161 

Schmidt G. D., Szkody P., Silvestri N. M., Cushing M. C, Liebert 
J., Smith P. S., 2005, ApJL, 630, L173 

SkrutskieM. F, etal., 2006, AJ, 131, 1163 

Steele P. R., Burleigh M. R., Farihi J., Giinsicke B. T., Jameson 
R. F, Dobbie P. D., Barstow M. A., 2009, A&A, 500, 1207 

Su K. Y. L., et al., 2007, ApJL, 657, L41 

Tremblay P.-E., Bergeron P., 2007, ApJ, 657, 1013 

Tremblay P.-E., Bergeron P., Gianninas A., 201 1, ApJ, 730, 128 

Vanlandingham K. M., Schmidt G. D., Eisenstein D. J., Harris 
H. C, Anderson S. F, Hall P. B., Liebert J., Schneider D. P., 
Silvestri N. M., Stinson G. S., Wolfe M. A., 2005, AJ, 130, 734 

Vennes S., Smith R. J., Boyle B. J., Croom S. M., Kawka A., 
Shanks T„ Miller L., Loaring N., 2002, MNRAS, 335, 673 

von Hippel T., Kuchner M. J., Kilic M., Mullally F, Reach W. T., 

2007, ApJ, 662, 544 

Vrba F. J., Henden A. A., Luginbuhl C. B., Guetter H. H, Munn 
J. A., Canzian B., Burgasser A. J., Kirkpatrick J. D., Fan X., 
Geballe T. R., Golimowski D. A., Knapp G. R., Leggett S. K., 
Schneider D. P., Brinkmann J., 2004, AJ, 127, 2948 

Wachter S., Hoard D. W., Hansen K. H, Wilcox R. E., Taylor 
H. M., Finkelstein S. L., 2003, ApJ, 586, 1356 

Warren S. J., et al., 2007a, ArXiv Astrophysics e-prints 

Warren S. J., et al., 2007b, MNRAS, 375, 213 

Wellhouse J. W., Hoard D. W., Howell S. B., Wachter S., Esin 
A. A., 2005, PASP, 117, 1378 

York D. G, et al., 2000, AJ, 120, 1579 

Zuckerman B., Becklin E. E., 1987, ApJL, 319, L99 

Zuckerman B., Koester D., Reid I. N., Hiinsch M., 2003, ApJ, 596, 
477 



APPENDIX A: WHITE DWARFS MODEL SPECTRA AND 
PHOTOMETRY 
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Figure Al. SDSS J003902.47-003000.3 model spectrum (black solid) 
with SDSS ugriz' (squares), UKIDSS YJHK (circles) and NTT JHK 
photometry (crosses). Also shown are the SDSS spectrum (light grey) and 
a composite WD+LO spectrum (dark grey). 



Figure A4. SDSS J013532.98+144555.8 model spectrum (black solid) 
with SDSS griz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+L5 
spectrum (dark grey). 
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Figure A2. SDSS J003923.04+003534.7 model spectrum (black solid) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+M9 
spectrum (dark grey). 



Figure AS. SDSS J075327.05+244705.5 model spectrum (black solid) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown is the SDSS spectrum (light grey) 
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Figure A3. SDSS J01 2032.27— 00 1351.1 model spectrum (black solid) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown is the SDSS spectrum (light grey) and an M7 spectrum scaled 
to the K-band magnitude of the putative secondary. 
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Figure A6. SDSS J085956.47+082607.5 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+M6 
spectrum (dark grey). 
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Figure A7. SDSS J090759.59+053649.7 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+L4 
spectrum (dark grey). 



Figure A10. SDSS J100259. 88+093950.0 model spectrum (solid black) 
with SDSS ugriz 1 (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+LO 
dwarf spectrum (dark grey). 
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Figure A8. SDSS J092648.48+102828.8 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+M8 
spectrum, (dark grey). 



Figure All. SDSS J 1 1 642.93+0443 17.7 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+M3 dwarf 
spectrum (dark grey). 




Figure A9. SDSS J09382 1 .34+342035.6 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+L5 
spectrum, (dark grey). 



Figure A12. SDSS J103448.92+005201.4 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+L5 dwarf 
spectrum (dark grey). 
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Figure A13. SDSS J103736.75+013912.2 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+M7 dwarf 
spectrum (dark grey). 



Figure A16. WD1 155-011 model spectrum (solid black) with SDSS 
ugriz' (squares) and UKIDSS YJHK photometry (circles). Also shown 
is a composite WD+M7 dwarf spectrum (dark grey). 
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Figure A14. SDSS Jl 10826.47+092721.5 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+M5 dwarf 
spectrum (dark grey). 
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Figure A17. SDSS J122625.93+002700.2 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+M6 dwarf 
spectrum (dark grey). 
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Figure A15. SDSS Jl 13416.09+055227.2 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+L3 dwarf 
spectrum (dark grey). 
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Figure A18. WD 1248+161 model spectrum (solid black) with SDSS 
ugriz' (squares) and UKIDSS JHK (circles) photometry. Also shown are 
the SDSS spectrum (light grey) and a composite WD+M8 dwarf spectrum 
(dark grey). The single J-band observation was classified as noise in the 
DR8 database and should be treated with caution. 
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Figure A19. WD1317+021 model spectrum (solid black) with SDSS 
ugriz' (squares) and UKIDSS YJHK photometry (circles). Also shown 
is a composite WD+L1 dwarf spectrum (dark grey). 



Figure A22. SDSS J134154.29+005600.8 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+M7 
dwarf spectrum (dark grey). 
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Figure A20. WD 1318+005 model spectrum (solid black) with SDSS 

ugriz' (squares), UKIDSS YJHK (circles) and AAO IRIS2 if -band pho- Figure A23. SDSS J141448.25+021257.7 model spectrum (solid black) 

tometry (cross). Also shown are the SDSS spectrum (light grey) and a com- with SDSS ugriz' (squares) and UKIDSS YJHK photometry (cir- 

bined WD + 600 K blackbody model (dark grey). cles). Also shown are the SDSS spectrum (light grey) and a combined 

WD+700K blackbody model (dashed grey). 
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Figure A21. SDSS J132925. 21 + 123025.4 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+M7 
dwarf spectrum (dark grey). 




1.5 

Wavelength (micron) 



Figure A24. SDSS J154431.47+060104.3 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+T3 
dwarf spectrum (dark grey). 
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Figure A25. SDSS J155720.77+091624.7 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey), a combined WD+L4 dwarf 
model (dark grey), and a combined WD+700 K blackbody model (dashed 
grey). 
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Figure A28. PHL 5038 model spectrum (solid black) with SDSS ugriz' 
(squares) and UKIDSS YJHK photometry (circles). Also shown are the 
SDSS spectrum (light grey) and a composite WD+L8 dwarf spectrum (dark 
grey). 
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Figure A26. SDSS J162514.88+302610.8 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+M5 
dwarf spectrum (dark grey). 
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Figure A27. SDSS J220841. 63-0005 14.5 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+L1 
dwarf spectrum (dark grey). 
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Figure A29. SDSS J22255 1 .65+001637.7 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS YJHK photometry (circles). 
Also shown are the SDSS spectrum (light grey) and a composite WD+L7.5 
dwarf spectrum (dark grey). 
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Figure A30. SDSS J233345.97-000843.0 model spectrum (solid black) 
with SDSS ugriz' (squares) and UKIDSS HK photometry (circles). Also 
shown are the SDSS spectrum (light grey) and a composite WD+M7 dwarf 
spectrum (dark grey). 



